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ABSTRACT 

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular 
clouds using data obtained with the PACS instrument onboard the Herschel Space Observatory and 
the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey 
(HOPS). A total of 55 new protostar candidates are detected at 70 jam and 160 /im that are either 
too faint (m24 > 7 mag) to be reliably classified as protostars or undetected in the Spitzer /MIPS 
24 /im band. We find that the 11 reddest protostar candidates with log AFx70/AFx24 > 1.65 are 
free of contamination and can thus be reliably explained as protostars. The remaining 44 sources 
have less extreme 70/24 colors, fainter 70 /im fiuxes, and higher levels of contamination. Taking the 
previously known sample of Spitzer protostars and the new sample together, we find 18 sources that 
have log AFx70/AFx24 > 1.65; we name these sources "PACS Bright Red sources", or PBRs. Our 
analysis reveals that the PBRs sample is composed of Class like sources characterized by very red 
SEDs (Tboi < 45 K) and large values of sub— millimeter fiuxes (Lgmm /Lboi > 0.6%). Modified black- 
body fits to the SEDs provide lower limits to the envelope masses of 0.2 Mq to 2 Mq and luminosities 
of 0.7 L0 to 10 L0. Based on these properties, and a comparison of the SEDs with radiative transfer 
models of protostars, we conclude that the PBRs are most likely extreme Class objects distinguished 
by higher than typical envelope densities and hence, high mass infall rates. 
Subject headings: ISM: Clouds, Stars: Formation 



1. INTRODUCTION 

The onset of the star formation process is broadly char- 
acterized by a dense collapsing cloud envelope surround- 
ing the nascent protostar. The dense cloud or protostel- 
lar envelope is opaque to radiation shortward of about 
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~ 10 jam and most of the radiation from these sources is 
reprocessed and emitted in the far-infrared (FIR). Fur- 
thermore, bipolar outflows from the protostar and disk 
carve out envelope cavities that enable a fraction of the 
protostellar luminosity to escape in the form of scattered 
light emission, predominantly at wavelengths shortward 
of ~ 10 jam. 

The earl iest phase of pro tostellar evolution, the Class 
phase ( An dre et al.l 119931 ), is thought to be short com- 
pared to the Class I phase (iLadal 11987). wit h combined 
Class and Class I lifetimes of ^0.5 Myr ( Evans et al.l 
120091 ): these estimates assume a constant star forma- 
tion rate and a typical age for Class II objects (pre- 
main sequence stars with disks) of 2 Myr. At the onset 
of collapse and immediately before the Class phase, 
protostars may go through a brief first hydrostatic core 
(FHSC) phase where the forming protostellar object be- 
comes opaque to its own radiation for the first time 
(|Larsonl Il969l ). The FHSC are expected to be very 
low-luminosity and deeply embedded. A population of 
very low luminosity protostars (Ve LLOs) were also re- 
cently identified by Spitzer (e.g., Dunham et al. 20081: 
iBourke et al. 2006), defined to have model-estimated in- 
ternal source luminosities of less than 0.1 Lq. VeLLOs, 
however, appear more evolved than FHSCs with features 
consistent with Class and I protostars. Furthermore, 
while several FHS C candidates have been id entified re- 
cently fe.g.. see [Enoch et al.l 120101: IChen et al. 20Tol: 
Pineda et al] 120111 : iPezzuto et al.l 1201 2[ ). it has proven 
observationally difficult to distinguish such sources from 
the young Class protostellar phase. It is therefore cur- 
rently difficult to identify the very earliest phases of the 
formation of a protostar. 
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Before the launch of Spitzer and the advent of ex- 
tremely sensitive mid-infrared surveys, conventional wis- 
dom held that a Class protostar should not be de- 
tectable at wavelengt hs shortward of 10 /ini due to 
the envelope opacity (Willia ms fc Ciezal l2Qllf ). Out- 
flows, however, can carve out cavities in the protostel- 
lar envelopes at a ver y early age and are e xpected to 
widen with evolution (jArce fc SargentI l2QQ6f ). Indeed, 
recent simulations have shown that even the extremely 
young FHSC sources may be capable of driving outflows 
(jCommercon et al. 2012; Price et al. 2012). Regardless 
of evolutionary state, the outflow cavities enable near- 
to mid-infrared light from the protostar and disk to es- 
cape and scatter off dust grains in the cavity or on the 
cavity walls. T his phenomenon has been well-k nown for 
Class I sources (iKenvon et al. II 19931 : [Padgett et al. 1999), 
but Class protostars were only well-detected in the 
mid-infrared wi th Spitzer (Noriega-Crespo et al.i i200^ 
iJg^rgensen et all 120071 : iStutz et al.ll2008D . The scattered 
light from Class protostars is often brightest at wave- 
lengths ~3.6 /im or 4.5 fim due to the dens e envelope 
obscurati on at shorter w avelengths (e.g.. iWhitnev et al.l 
^003a; To bin et al.ll2007[ ). 

The combination of results from recent space- 
based and ground-based surveys have resulted in 
well-sampled spectral energy distributions from the 
near-infrared to the (sub) millimeter fo r large sam- 
ples of protostellar objects (e.g., Hatche lT et al.l [20071: 
Enoch et al^ '2009'; 'Launhardt et aL| 120101 : [Fischer et all 
20101 : [Lminhardt et al. 2012). These SEDs are domi- 
nated by scattered light between ~ 1 /im and 10 /im, opti- 
cally thick thermal dust emission from ~ 10 to ~ 160 /im, 
and optically thin dust emission at wavelengths longward 
of ~ 160 /im. Radiative transfer models of protostellar 
collapse have become increasingly important to interpret 
these data since these can account for the varying tem- 
perature and density p rofiles in the envelope s surround- 
ing the protostar (e.g.. IWhitnev et alll2003al lb[). 

The large number of free model parameters — such 
as the combinatio n of outflow cavities, rotationally flat- 
tened env elopes ( Ulrichl ll976: Cassen fc Moosman '1981| 
iTerebev et al. 198^ . and varying viewing angles — can 
however make th e best-fit SEP model p arameters highly 
degenerate (e.g.. IWhitnev et aH l2003al ). For example, 
sources viewed at nearly edge-on orientations can be sub- 
stantially more obscured than sources at the same evolu- 
tionary state viewed from a less extreme vantage point. 
Thus, standard diagnostics such as bolometric tempera- 
ture or mid-infrared spectral index can yield vastly dif- 
ferent results depe nding on the source inclination (e.g., 
[Dunham et al.|[2Q1 0). Radiative transfer models can help 
break some of these degeneracies but ambiguities can re- 
main as to whether a source has a very dense envelope 
or if it is simply viewed edge-on. 

While much has been learned about the Class phase 
from observations and modeling, there are relatively few 
Class object s present in ne arby star-forming clouds and 
globules (Evans e t al.[[2009l ) compared to the numbers of 
Class I and Class II sources. One of the principal goals 
of recent star formation surveys has bee n to understand 
the ev olution of protostellar sources. [Young fc Evans[ 
(|2005[ ) generated models for the smooth luminosity evo- 
lution of protostellar objects that will become 0.3 
1 Mq and 3 Mq stars. These models, however, over 



predict the luminosities of most protostars located in 
nearby star-forming regions, a f act that is taken as 
evidence for episodic accretio n ([Kenyon fc Hartmand 
1990; Evans et al. 2009; Dunha m et al.[[2010[ ). However. 
lOffner fc McKed ([2011) show that the observed luminos- 
ity functions of protostars can be explained through a 
dependence of the mass accretion rate on the instanta- 
neous and final mass of the protostar. The low resolu- 
tions and sensitivities of previous FIR instrumentation 
have made the detection of protostars in more distant 
and richer star-forming regions difficult and subject to 
substantial confusion. Thus, studies of protostellar evo- 
lution have been limited to combining all known Class 
proto stars from the nea r by regions into a s ingle analysis 
(e.g., [Mvers et all[1998[ : [Evans et al1[2009[ ) to achieve a 
more robust sample size. 

The advent of the Herschel Space Observatory 
(i Pilbratt et al.[ [2QTq[) has tremendously improved reso- 
lution and sensitivity to FIR radiation, where protostars 
emit the bulk of their energy. These improvements en- 
able the study of protostellar populations to be extended 
to more distant, richer regions of star formation that have 
more statistically significant samples of prot ostars in the 
Class and I phases (e.g., Ragan et al.[[2012[ ). The Her- 
schel Orion Protostar Survey (HOPS) is a i/"ergc/ie/ Open - 
Time Key Programme (OTKP) (e.g . jStanke et"al][201C ; 
"Fischer et alj \201il [Ali et al.[ [2010t [Manoj et all [2012 ) 
targeting ~ 3 00 of the Spitzer ident ified Orion protostars 
with PACS (Poglitsch et aU [20101 ) 70 jam and 160 /im 
photometry and P ACS spectroscopy (53 /im to 200 /im; 
[Manoj et al.[[2012[ ) for a subset of 30 protostars. Orion 
is th e richest star-forming region within 500 pc of the 
Sun ([Megeath et al.[ [2012|). a t a dista nce of ^420 pc 
(average value: Mente n et al.[ [2007c iHirota et al.[ [2007[ : 
Sandstrom et al. 2007). The large sample of protostars 
in Orion enables studies of protostellar evolution to be 
carried out for a single star forming complex where all 
protostars lie at nearly the same distance with a statisti- 
cally significant sample, comparable to or larger than all 
the nearby regions combined. The large sa mple may also 
enable short timescale phenomena (e.g., iFischer et al.[ 
[2012) to be detected such as brief periods of high en- 
velope infall rate in the earliest phases of star formation. 
Even with the increased numbers of protostars in Orion, 
however, it is unclear if we would expect to detect FHSCs 
given the faintness of these sources and short lifetimes of 
less than 10 kyr (Commerco n et aP 2012). 

The PACS imaging of the HOPS program has the po- 
tential to identify protostars that were not detected by 
Spitzer due to a combination of opacity of the envelope 
and/or confusion with nearby sources. Indeed, the PACS 
70 /im band of Herschel is ideal for detecting such pro- 
tostars, with the highest angular resolution, limiting the 
blending of sources. Also, the lower opacity relative to 
MIPS 24 /im allows the reprocessed warm inner enve- 
lope radiation to escape. Finally, and most importantly, 
a 70 /im point source is strong evidence for an embedded 
protostar because external heating cannot raise temper- 
atures high enough to emit at this wavelength. Thus, 
some cores in sub-millimeter surveys that were previ- 
ously identified as starless may in reality be protostellar. 

Using Herschel^ we have serendipitously identified a 
sample of 70 jam point sources that were not identified 
in the previous Spitzer protostar sample (Megeath et al.[ 
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l2Q12f ). Furthermore, we have identified a subset of these 
that have the reddest 70 fim to 24 jam colors of ah proto- 
stars in the combined Orion sample. These sources may 
have the densest envelopes and are possibly the youngest 
detected Orion protostars and we name them "PACS 
Bright Red sources" or PBRs. 

We will describe the methodology of identifying these 
sources and classify them as either being protostellar, 
extragalactic contamination, or spurious detections co- 
incident with extended emission. We will discuss the ob- 
servations and data reduction in § 2, the source finding 
and classification methods in § 3, the observed properties 
of the new sources in § 4, the PBRs in § 5, the compar- 
ison of the cold PBRs to models in § 6, some relevant 
model degeneracies in § 7, and, finally, our results in 
§ 8. Throughout this work, all positions are given in the 
J2000 system. 

2. OBSERVATIONS, DATA REDUCTION, AND 
PHOTOMETRY 

In this work, we present Herschel scan map observa- 
tions of a sub-set of the HOPS fields containing can- 
didate protostars. In addition, we present a subset of 
our APEX LABOCA and SABOCA observations of these 
fields. A summary of the HOPS Herschel PACS survey 
observations is presented in Tables [1] and O Here we dis- 
cuss the observations, data processing, and photometry 
extraction. 

2.1. Herschel PACS 

The PACS data were acquired simultaneously at 70 fim 
and 160 /im over 5' x 5' or 8' x 8' field sizes. The field sizes 
and centers were chosen to maximize observing efficiency 
by allowing each field to include as many of the target 
Spitzer-selected protostars (Megeath et al. 2012) as pos- 
sible while minimizing redundant coverage. The obser- 
vations were acquired at medium scan speed (20^^ sec~^), 
and are composed of two orthogonal scans with homoge- 
neous coverage. 

The PACS data were reduced using the Herschel Inter- 
active Processing Environment (HIPE) version 8.0 build 
248 and 9.0 build 215. We used a custom built pipeline to 
process data from their raw form (the so-called Level 
data) to fully calibrated time lines (Level 1) just prior to 
the map-making step. Our pi peline uses the same pro- 
cessing steps as described by Poglitsc h et al.l jJoTo) but 
also include the following additions and modifications. 
First, we used a spatial redundancy based algorithm to 
identify and mask cosmic ray hits. Second, we mitigated 
instrument cross-talk artifacts by masking (fiagging as 
unusable) detector array columns affected by cross-talk 
noise. This technique is effective but at the expense of 
loss of signal from the affected detector array columns. 
Third, we used the "FM6" version of the instrument re- 
sponsivity, which has a direct bearing on the absolute 
calibration of the final mosaics. 

The Level 1 data were processed with "Scanamorphos" 

ousse 1(2012) version 14.0. The final maps were pro- 
duced using the galactic option and included the turn- 
around (non-zero acceleration) data. The final map 
pixel scales used in this work are 1.0'Ypix at 70 /im and 
2.0'7pix at 160 /im. 

The photometry was performed in the following fash- 
ion. We first derived customized aperture corrections 



to the 70 /im and 160 /im data using the Herschel Sci- 
ence Center (HSC) provided observations of Vesta (to 
be discussed in more detail in Fischer et al., in prepa- 
ration). At 70 /im, we used radii of sizes 9.6'^ 9.6'^ 
and 19.2^', for the aperture and sky annuli respectively. 
For these parameters we derived an aperture correc- 
tion of 0.7331, where the measured fiux in the aperture 
is divided by this correction to obtain a total point- 
source fiux. At 160 /im, we used aperture radii of 12. 8'^ 
12. 8'^ and 25. 6'^ for the aperture and sky annuli respec- 
tively. Similarly, we derived an aperture correction of 
0.6602. The encircled energy fractions provided by the 
photApertureCorrectionPointSource task in HIPE do 
not account for the effect of applying an inner sky an- 
nulus that is close to the size of the source aperture and 
small compare to the PSF. Our apertures therefore ac- 
count for 3-4% of the source ffux that is removed. Fur- 
thermore, our adopted aperture sizes are smaller than the 
PACS instrument team recommendation but were chosen 
to minimize contribution from nebulosity (extended, non 
point-like emission) often surrounding the protostars in 
Orion. Given the complex structure in the images and 
at times crowded fields, our aperture photometry may 
suffer from blending and contamination. The photomet- 
ric errors include a 10% systematic error floor added in 
quadrature to the standard photometric uncertainties. 
These errors represent systematic uncertainties in our 
photometry and aperture correction, as well as the over- 
all calibration uncertainty of PACS. We note that the 
reported HSC point-source calibration uncertainties for 
PACS are ~ 3% at 70 /im and ~ 5% at 160 /im, and were 
derived from isolated photometric standards. Therefore, 
our flnal uncertainties are conservative. 

We inclu de the 100 /im Gould Belt S urvey 
(GBS; e.g., lAndre et all [20Tol : iKonvves et all [20T0I : 
iMen'shchikov et all [2QTqI . see also Schneider et al., in 
preparation, for Orion B, and Roy et al., in preparation, 
and Polychroni et al., in preparation, for Orion A) data 
of Orion in this work for the PBRs analysis. Given the 
sparsely covered SEDs of our sources, these data pro- 
vide important information regarding the shape of the 
thermal SED of cold envelope sources. These data were 
acquired using the medium scan speed (20''sec~^) and 
cover an area much larger than the HOPS flelds. These 
data were processed in a similar way to the HOPS pro- 
cessing described above. Following the above 70 /im and 
160 /im analysis, we used aperture radii of sizes 9.6'^ 
9.6'^ and 19. 2'^, for the aperture and sky annuli respec- 
tively. For these parameters, we derived an aperture cor- 
rection of 0.6944. As with the 70 /im and 160 /im data, 
we also assume a conservative 10% systematic error floor. 

2.2. APEX SABOCA and LABOCA 

We obtained sub-millimeter (smm) continuum maps 
using the LABOCA and SABOCA bolometer array s on 
the APEX telescope. LABOCA (Siringo et aiJ|2009j) is a 
~250 bolometer array operating at 870 /im, with a spa- 
tial resolution of ~ 19'^ at FWHM. We used a combina- 
tion of spiral and straight on-the-fly scans to recover ex- 
tended emissio n. Data reduc tion was done with the BOA 
software (^Schuller et al .11201 2) following standard proce- 
dures, including it erative source modeling. SABOCA 
(tSiringo et aD[2QlQD is a 37 bolometer array operating at 
350 /im, with a resolution of -7.3'' FWHM. The observ- 
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ing and data reduction procedures were similar to those 
used for LABOCA. For both cameras, observations were 
carried out between 2009 November and 2012 June, and 
are still ongoing to complete our Submillimeter Orion 
Survey. Conditions were generally fair over the course of 
our observing campaign. The observations will be sum- 
marized in more detail by Stanke et al. 2013, in prepa- 
ration. 

The beam sizes of the final reduced maps are 7.34'^ 
and 19.0'' FWHM for the SABOCA and LABOCA ob- 
servations, respectively. The photometry was extracted 
in the same way for both wavelengths. When possible, if 
there was a strong source detection, we re-centered using 
the 70 jam catalog source coordinates. Given the con- 
tributions of flux due to surrounding cold material, such 
as fllaments and other extended envelope structure, it is 
likely that a single photometric measure can suffer from 
large systematic effects. We have measured source fluxes 
in three ways: 

1. We measured the source peak flux per beam. 

2. We measured source flux over an aperture with radius 
equal to the FWHM at the corresponding wavelength (r 
= 7.34'' and 19.0" at 350 and 870 /im, respectively), us- 
ing a sky annulus with inner and outer radii equal to 
[1.5,2.0]xFWHM, corresponding to 11.0" and 14.7" at 
350 /im, and 28.5" and 38.0" at 870 jam. 

3. We measured the flux over the same aperture size as 
the previous method without any sky subtraction. In the 
case that a source was not strongly detected and we were 
not able to re-center, the 70 /im catalog source coordi- 
nates was used, along with method 3, and the photomet- 
ric point was flagged as an upper limit. By re-centering 
whenever possible, we accounted for possible pointing 
offsets between data-sets, which can be signiflcant. The 
calibration error dominated the error budget for well de- 
tected sources; we therefore adopted a flux error equal 
to 20% and 40% of the measured flux for LABOCA and 
SABOCA respectively. The photometric fluxes are pre- 
sented in Table [5l 

2.3. Spitzer IRAC and MIPS 

The IRAC and MIPS imaging and photometry pre- 
sented here are taken from the 9 degree^ survey of 
the Orion A and B cloud obtained during the cryo- 
genic Spitzer mission. The data analysis, extraction 
of the IRAC 3.6 /im, 4.5 /im, 5.8 /im, and 8 /im and 
MIPS 24 /im photometry, and the compilation of a point 
source catalog containing the combined 2MASS, IRAC 
and M IPS photo metry are descr i bed in iMegeath et al.i 
(|2012f ): also see iKrvukova et all (j2012D for a detailed 
description of the MIPS 24 jam photometry. In total, 
298405 point sources were detected in at least one of the 
Spitzer bands, and 8021 sources were detected at 24 /im 
with uncertainties < 0.25 mag. The Spitzer images used 
in this work are taken from the mosaics generated from 
the Orion S urvey data using Clus ter Grinder for the 
IRAC data (|Gutermuth et al.! l2009h and the MIPS in- 
strument team's Data Analysis Tool for the 24 /im data 
(Gordon et al. 2005). The MIPS data are saturated to- 
ward the Orion Nebula and parts in the NGC2024 region; 
we exclude these saturated regions from our analysis. 

The identiflcation of protostars with the Spitzer data 
was based primarily on the presence of a flat or rising 
spectral energy distributions between 4.5 /im and 24 /im 



(|Krvukova et al.ll2012l : rMegeath et alll2012f ). In addition, 
Megeath et al. identified objects which have point source 
detections only at 24 jam but which also showed other 
indicators of protostellar nature such as the presence of 
jets in the IRAC bands. To minimize contamination from 
galaxies, the Spitzer-identifLed protostars were required 
to have 24 /im magnitudes brighter than 7th magnitude; 
fainter than 7 magnitudes, the number of background 
galaxies begins to domi nate over t he number of embed- 
ded sources (Krvukova et al.l [20121 ). Given the imposed 
24 /im magnitude threshold, the faintest and reddest 
protostars may not be inclu ded in the Spitzer sample. 
In total, the Megea th et al . ('2012') catalog contains 488 
protostars. Of these, 428 are classified as bona fide pro- 
tostars, 50 are faint candidate protostars, and 10 are red 
candidate protostars. The faint candidate protostars are 
sources with 24 /im magnitudes higher than 7.0. The red 
protostars are sources that are only detected as a point 
source at 24 jam and are thus not classifiable. Due to 
their location within high extinction regions and/or as- 
sociation with jets or compact scattered light nebulae in 
the IRAC bands, , they have been included in the cata- 
log. (Indeed, this last category was added to the Megeath 
et al. catalog after the Herschel data revealed that a 
relatively large number of such sources would likely be 
confirmed as protostars.) In addition, the Megeath et 
al. catalog identified 2992 objects as pre-main sequence 
stars with disks. In what follows, we use the protostar 
catalog of 488 Spitzer sources to catalog previously iden- 
tified sources in the HOPS images. 

In contrast to the full Megeath et al. catalog, the HOPS 
protostar sample represents a sub-set of that catalog. 
The HOPS sample is not a complete representation of 
the Megeath et al. catalog because the HOPS survey tar- 
geted a smaller area then the Spitzer survey and required 
an extrapolated 70 /im fiux that would be detectable 
with PACS. However, the HOPS catalog represents our 
best pie-Herschel knowledge of the protostellar content 
within the HOPS fields. The HOPS sample is therefore 
the best one against which to bench-mark and compare 
the new Herschel protostellar candidates. The HOPS 
sample constitutes ~ 300 previously identified Spitzer 
protostars all having Herschel PACS 70 /im detections. 
Of these, - 250 HOPS sources have PACS 160 /im detec- 
tions. This sample will be analyzed in detail in Fischer 
et al., in preparation. 

In contrast to the full Megeath et al. catalog, the 
HOPS protostar sample are those protostars specifically 
targeted by the HOPS program. The majority of the 
HOPS sample consists of S pit zer-identifi.ed protostars 
with 24 /im detections; hence, protostars in regions that 
are saturated in the 24 /im images of Orion, namely the 
brightest regions of the Orion Nebula and NGC2024, are 
not included. These protostars were also required to have 
a predicted 70 jam fiux > 20 mJy as extrapolated from 
their 3 jam to 24 /im SEDs. In addition, protostar candi- 
dates with only 24 /im detections were included if there 
was independent information of their protostellar nature. 
The HOPS sample represents the best pie-Herschel cat- 
alog of protostars that were expected to be detected with 
Herschel /F ACS and were not found in bright nebulous 
regions. There are ~ 300 protostars in the HOPS catalog 
which have been detected at 70 jam and ~ 250 protostars 
detected at 160 um. The uncertainty in there absolute 
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Fig. 1. — Frequency of 70 /im FWHM values for all sources de- 
tected in our HOPS images that have not been previously identi- 
fied in the Spitzer catalog. The black line shows our adopted 7.8^^ 
FWHM threshold, above which we reject all sources. Inset: Same 
as above for the subset of sources with both FWHM (70 /im) < 7.8^^ 
and 160 /im detections (see § 3). The dotted histogram indicates 
the 160 /im FWHM distribution. 



number of protostars is due to the ongoing process of 
eliminating contamination from the sample. 



3. 



IDENTIFICATION OF NEW CANDIDATE Herschel 
PROTOSTARS 



To find protostars which were not reliably identified 
with Spitzer^ we must first isolate a sample of sources 
that are detected in the PACS 70 /im band but are ei- 
ther fainter than 7.0 magnitudes or undetected in MIPS 
24 /im waveband. To identify such sources in each HOPS 
field we first g enerate a 70 /im sourc e catalog using the 
Phot Vis tool (IGutermuth et al.ll2008f ). The Phot Vis tool 
uses a sunken Gaussian filtering to extract sources that 
are of order the size of the Gaussian FWHM, an input 
parameter. We choose this parameter to be the size of 
the 70 /im PSF FWHM, or 5". PhotVis also requires 
a SNR threshold as input; we adopt a low value of 7 
to balance the recovery of as many candidate sources as 
possible while still rejecting noise spikes. 

Furthermore, we must reject unreliable sources near 
the edges of maps where the lower coverage causes el- 
evated noise levels. The Scanamorphos scan map im- 
age cubes include a weight map for each field. Within 
Scanamorphos, the weight map is computed over the 
same projection as the sky map, a nd is defined as one 
over the variance in the white noise (jRousselll20l"2h . Each 
weight map is then normalized by the average map value 
(fRoussel 2012). We find that for the HOPS data-set, 
weight map values of '^20 are confined to the outer higher 
noise edges of our scan maps. We therefore use the weight 
maps to reject edge sources from the catalog at this phase 
of the analysis. We accomplish this by requiring that the 
mean value of the weight map in a 9 x 9 pixel area cen- 
tered on the candidate source have a value of at least 20. 
For reference, all HOPS 70 /im scan maps have weight 
map values greater than 60 over most of the map areas. 

The resulting preliminary source catalog includes all 
sources in the 70 /im images, i.e., previously identified 
Spitzer sources, new candidate protostars, nebulosity, 
and other undesirable features and artifacts in the im- 



ages. We then cross-correlate this PACS 70 /im prelim- 
inary catalog with the existing Spitzer catalog to elimi- 
nate all previously identified protostars in each field that 
are brighter than the previously adopted 24 /im cutoff 
of 7 mag (Megeath et al. 2012). Therefore our sample 
includes by definition only sources that are faint or un- 
detected in the previous Spitzer catalog. 

To identify previous source detections, we require that 
a source be matched to within a positional offset of 8'^ 
when cross-correlated with the Spitzer catalog. This 
threshold is conservatively large compared to the Spitzer 
astrometry and is meant to encompass two main sources 
of astrometric error. First, it is possible that the absolute 
coordinates of a source may shift as function of wave- 
length (although this effect is expected to be relatively 
small) since different wavelengths may trace different ma- 
terial near the protostars. Second, the Herschel point- 
ing accuracy, which is of order ^ 2" (Icr), dominates 
the positional accuracy for most sources when compar- 
ing the Spitzer catalog source coordinates to the Herschel 
70 /im coordinates. To match coordinates robustly, we 
therefore adopt a conservatively large 8'^ threshold. We 
find that inspection of the matched sources by eye shows 
that this threshold works well and provides a low rate of 
mismatched or duplicate sources. 

Our final goal is to obtain a sample of previously 
unidentified and uncharacterized Herschel protostar can- 
didates. These sources should be characterized by a 
point-like appearance at 70 /im. Therefore, after reject- 
ing all Spitzer protostars as described above, we next 
apply a simple FWHM (or apparent size) filter to the 
remaining 70 /im sources. The distribution of 70 /im az- 
imuthally averaged FWHM values is shown in Figure [T] 
as the solid black histogram. We find a clear peak in the 
distribution at low FWHM values, indicating a popula- 
tion of point-like sources. Based on this distribution, we 
adopt a FWHM threshold of 7.8'^ meant to select 70 /im 
point-sources. We find 127 sources that fulfill the crite- 
ria listed here: 85 of these have 24 /im detections while 
42 do not. 

In a further step, we then require that all 70 /im sources 
also have a 160 /im detection and not upper limits. The 
70 /im FWHM distribution of this sub-set of sources is 
shown in the inset of Figure [T] along with the correspond- 
ing 160 /im FWHM distribution. Our final sample con- 
sists of 55 candidate Herschel protostars with both 70 /im 
and 160 /im detections. Of these, 34 have Spitzer 24 /im 
detections fainter than 7.0 magnitudes and 21 do not 
have any 24 /im detection. 

3.1. Spitzer non-detections 

A search for newly detected protostars using Herschel 
requires us to determine upper limits at 24 /im for those 
sources that are not detected by Spitzer. To determine 
these limits, we adapt the method developed by Megeath 
et al. to assess the spatially varying completeness of the 
Spitzer Orion Survey data. The completeness of the 
24 /im data depends strongly on the presence of neb- 
ulosity and point source crowding. To account for these 
factors, we measure the fiuctuations of the 24 /im signal 
in an annulus centered on the position of the Herschel 
point-source using the the root medi an square deviation , 
or RMEDSQ (see Equation 1 of Megeath etaD |2012D . 
We then use the results from the artificial star tests (see 
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Fig. 2. — IRAC color indices 0:3.6-4.5 and 0:5.8-8.0 for sources 
detected in at least one IRAC band. The Herschel protostar can- 
didate sample is shown as squares (sources with 24 fim detections) 
and circles (sources without 24 fim detections), while the HOPS 
protostar sample is shown as triangles. Thick bold points are 
sources with robust 870 /im detections. Candidate protostars with 
<^5. 8-8.0 > 3.0 or upper limits consistent with this threshold are 
flagged as possible extragalactic contamination and highlighted in 
red. Note that some candidate protostars have both 05.8-8.0 > 3.0 
and a robust 870 /im detection; these sources are not considered 
extragalactic contamination. 
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Fig. 3. — IRAC and PACS 160 /im color indices 03.6-4.5 and 
03.6— 160 for the candidate protostar sample. Sources previously 
characterized as extragalactic based on their 05.8-8.0 > 3 index (or 
limit) are indicated in red with a "G". Sources with 05.8—8.0 < 3 
are indicated in orange. The orange lines indicate the lower limits 
of 03.6— 160 for sources with no IRAC detections. Sources with 
03.6— 160 > 1-2 are considered highly probable protostars, while 
other sources are flagged as described in the text. 



Appendix of lMegeath et alJl2Q12f ) to determine the mag- 
nitude at which 90% of the point sources would be de- 
tected for the observed level of fluctuations. We convert 
this magnitude into a flux density to obtain 24 /im upper 
limits. 

Several of the identifled protostars show IRAC emis- 
sion but are not included in the Megeath et al. point 
source catalog because they are spatially extended. To 
obtain homogeneously extracted IRAC fluxes for the en- 
tire sample of sources, we measure fluxes using an aper- 



ture of 2 pixels, with a sky annulus of 2 to 6 pixels, cor- 
responding to aperture radii of 2AA" ^ 2AA" ^ and 7.33'^ 
respectively, with a pixel scale of 1.22'^ pixel" ^. We use 
the PACS 70 /im source coordinates as starting guesses, 
and attempt to re-center at each IRAC wavelength. If 
the re-centering fails, as for sources with no IRAC de- 
tections, we take the integrated flux in that aperture at 
the original PACS 70 /im source coordinate to be the 
upper limit. The aperture corrections and photometric 
zero-points are those given by Kryukova et al. 

3.2. Contamination in the sample 

Galaxies often exhibit infrared colors similar to those of 
young stellar objects (YSOs) due to the presence o f dust 
and hydrocarbons in the galaxies (St ern et"ani2005[ ). Ex- 
tensive work has been done towards characterizing the 
extragalactic "contamination" in Spitzer surveys of star- 
forming regions and mitigating it through photomet- 
ric crit eria designed to separate galaxies from bona fid e 
YSOs (jGutermuth et al.l[2009l . l2008l : lHarvev et al.ll2007D . 
These authors show that star-forming galaxies can be 
distinguished from YSOs by the galaxies' stellar-like 
emission in the IRAC 3.6 jam and 4.5 /im bands and their 
bright polycyclic aromatic hydrocarbon (PAH) emission 
in the IRAC 5.8 /im and 8.0 /am bands CG utermu th et al.l 
l2009l . l200l : IWin^ et al]l2007l : IStern et al. 2005). How- 
ever, we note that some Active Galactic Nuclei (AGN) 
dominated galaxies may not exhibit PAH emission; 
therefore an analysis based only on the IRAC colors 
may not capture all p ossible sources of contamination 
(jRobitaille et al.ll2008l ). 

To analyze the IRAC colors of our sample, we define 
a = dlog{XFx) / dlog{X) . In Figure O we plot a^,8-8 
vs. 0^3.6-4.5 for the sample of new Herschel sources with 
coverage in all four of the Spitzer /IRAC bands and de- 
tections in at least one band, compared to the HOPS 
protostar sample. This color index is relatively insensi- 
tive to reddening since the extinction in the 5.8 /im and 
8 /im band s of IRAC are very similar (e.g. jFlahertv et all 
2007; Guter muth et al. l 2008). Figure [2] shows a cluster of 
sources with high values of 0^5. s-s (i-e., 0^5.8-8 ^ 3; solid 
horizontal line) yet 0^3.6-4.5 values of an SED that is de- 
clining or fiat with increasing wavelength. These sources 
show the characteristics of star-forming galaxies with 
bright PAH emission (resulting in high values of 0^5. 8-8) 
but values of 0^3.6-4.5 that are dominated by starlight. In 
our adopted scheme, 0^5.8-8 ^ 3 corresponds to a color 
of [5.8] - [8] > 2.17; this thre shold is higjher than the 
[5.8] - [8] > 1 threshold used bv lGutermuth et al.l (|2QQ8[ ) 
to isolate galaxies and thus ensures that most protostellar 
candidates will be less likely to be miss-identified extra- 
galactic sources (Allen et al. 2004; Megeath et al. 20dl). 
We identify the cluster of sources with 0^5.8-8 > 3 and 
<^3.6-4.5 < 0.5 as likely extragalactic contamination. We 
note that nebular contamination of the photometry can 
cause PAH-like 0^5. 8- 8 values, and thus may cause us to 
over-estimate the extragalactic contamination. Of the 
55 sources identified here as protostellar candidates we 
flag 23 as possible extragalactic contamination based on 
this criteria. However, other sources of contamination, 
such as AGN lacking PAH emission, may remain in our 
sample. 

Inspection of the SEDs of the remaining 32 sources 
show a range of SED slopes and shapes. It is possi- 
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Fig. 4. — Example SEDs for each of the 5 categories described in § 3.2: green SEDs (flag = 1) are considered to be rehable protostar 
candidates, yellow SEDs (flag = 2) are less likely to be protostellar, red SEDs (flag = 3) are considered extragalactic contamination based 
on the 8 /im PAH feature, the purple SED (flag = 4) is an example of one of 3 sources that may be explained as transition disks, and, 
flnally, the black SED (flag = 5) is unclassiflable due to the lack of IRAC coverage. Errors are smaller than the size of the points, except 
for the 350 /xm and 870 /im points; upper limits are indicated with triangles. 



ble that such sources may also be extragalactic contam- 
ination by AGN, for example. To address this, we fur- 
ther refine the analysis presented above by analyzing the 
3.6 fim to 160 /im SED shapes with the spectral index 
<^3. 6-160- As illustrated in Figure [3l the sources fiagged 
as extragalactic based on the 0^5. s-s index (red points) 
generally have 0^3.6-160 ^ 1-2. We therefore calibrate 
the 0^3.6- 160 relative the reliable extragalactic candidates 
with robust IRAC detections and expect that extragalac- 
tic sources will have 0^3.6-160 ^ 1-2. Using this crite- 
rion, we refine our source classification as follows. All 
sources with 0^3.6-160 > 1-2 (and 0^5. 8-8 < 3 when IRAC 
detections exist) are fiagged as high probability proto- 
stars (fiag = 1 in Table [3j). Sources having values of 
0-5 ^ <^3. 6-160 ^ 1-2 but that originally classified as 
candidate protostars based on a low value of 0^5. 8-8 are 
fiagged as less likely to be of a protostellar nature (fiag 
= 2 in Table [3]). Furthermore, by definition, sources 
originally classified as extragalactic based on their PAH 
signature at 8 /am remain classified as such (fiag = 3 in 
Table [3]). Sources with 0^3.6-160 < 0.3 and 0^3.6-4.5 ~ —3 
are fiagged as "other" (fiag = 4 in Table 3) since their 
SEDs are consistent with a stellar photosphere at shorter 
wavelengths. Finally, one source has no IRAC coverage 
and therefore is fiagged with a value of 5. In Figure 2] we 
show example SEDs of each category. 

Only one source (313006) originally fiagged as ex- 
tragalactic based on its 0^5.8-8 limit (non-detection at 
5.8 /im and 8 /im) was revised to a highly probable pro- 
tostar (see Figure[3]and top left panel of FigureH]). In ad- 
dition, as we note above, we find three sources with SEDs 
that we label "other" (fiag value = 4) which are inconsis- 
tent with the categories described above. Source 069001 
(see Fi gure [H top r ight p anel) was previously character- 
ized bv lFang et al.l (|2009f ) as a K7 star with a debris disk, 
with a very poorly constrained age of ~ 0.061q;q3 Myr. 
These authors only include data up to 24 /im. The SED 



we observe with Herschel may be consistent with a tran- 
sition disk but not a debris disk. The remaining two 
sources in this category have similar SEDs as that of 
069001; while a transition disk explanation for all three 
sources may appear likely depending on the age of the 
sources, we cannot currently rule out other possibilities. 
Nevertheless, all of these sources have SEDs consistent 
with a stellar photosphere in the IRAC bands, and hence 
these are likely to be fully formed stars surrounded by 
circumstellar dust. 

Interestingly, we find that the most reliable SED clas- 
sification criterion by far is that of sources that have 
neither IRAC nor MIPS 24 /im detections. Of these, we 
find that all 6 sources have strong sub-millimeter detec- 
tions and reside in dense and filamentary environments. 
This finding points to the critical importance of obtain- 
ing high resolution sub-millimeter data to constrain the 
properties of such sources. In the following text, we in- 
clude all 55 Herschel-detected sources in our analysis 
and figures. 

4. Herschel PROTOSTAR CANDIDATES 

We present the Herschel protostar candidate catalog in 
Table [3l Here we include the PACS 70 jam coordinates 
and fiux measurements at 24 /im, 70 /im, and 160 /im. 
We indicate which sources are fiagged as reliable pro- 
tostellar candidates and which are likely contamination, 
based on the results from the previous section. We also 
indicate if the sources have a robust 870 jam detection. 
Furthermore, we present the values of Lboi and Tboi and 
their corresponding estimated statistical errors (see dis- 
cussion in § 5.1). In Figure [51 we show the 70 /im fiux 
distributions for the sample compared to the distribution 
of HOPS protostars. The majority of the new candidate 
protostars have 70 jam fiuxes that are lower than the 
previously identified Spitzer HOPS sample; this is not 
surprising since the new candidate protostar sample is 
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Fig. 5. — Distributions of 70 iim fluxes for the HOPS protostars 
(dashed histogram), new candidate protostars (sohd histogram), 
and the 18 reddest sources drawn from both combined samples 
(dot-dashed histogram). The corresponding grey hues indicate the 
median flux values (vertical lines) and the 68% interval (horizontal 
lines). 



selected to be faint or undetected at 24 jam. Further- 
more, the peak at low 70 jam flux values is dominated by 
extragalactic contamination, as discussed above. 

In Figure El we show the MIPS 24 /im, PACS 70 jam, 
and PACS 160 /im colors of the new Herschel sources 
compared to the HOPS sample of 70 /im detected pro- 
tostars. The top panel shows the 70 /im flux vs. the 
log (AFa70/AFa24) color (henceforth 70/24 color), while 
the bottom panel shows the log (AFxl60/AFx70) color 
(henceforth 160/70 color) vs. the 70/24 color, for our 
sample of new protostar candidates compared to the col- 
ors of the Spitzer-identifLed HOPS sample. The Spitzer 
24 /im 7 magnitude limit, imposed on the HOPS sam- 
ple for a reliable protostellar identification, is apparent 
in the top panel as the diagonal line approximately sep- 
arating the new protostellar candidates at fainter 70 /im 
fluxes and redder 70/24 colors from the population of 
Spitzer-identifLed HOPS sources. 

For comparison, in the top panel of Figure [6l we 
also show the fluxes and colors of presumably typi- 
cal and well-studied Class sources: VLA1632-243 (J. 
Green and DIGIT team, private communi cation, 2012, 
and Green e t al., in pre p. ), IRAS16293 (lEvans et"aLl 
l200a . B335 (IStutz et al.l 120081 : iLaunhardt et al.ll2012D 
CB68 (iLaunhardt et al.ll2012f ). and CB244 (jStutz et al.l 
l2010l : [Launhardt et al.ll2012f ). Furthermore , we also show 
the c olors of various VeLLOs: L673-7 ( i Dunh am et al. 
[20081), IRAM04191 (iDunham et al.l I2006D . and CB130 
(|Launhardt et al.ll2012 l. We find that the observed col- 
ors of our sample of candidate protostars appear consis- 
tent with the colors of more near-by Class and VeLLO 
sources but not with FHSC candidate colors proposed 
in the literature (e.g.. iCommercon et al.ll2012f ). We find 
that the majority of these previously known Class and 
VeLLO sources do not appear as red in their 70/24 colors 
as the reddest sources in our sample. The only exceptions 
to this trend are IRAS16293 and VLA1632-243, perhaps 
representing an extrema in the 70/24 color distribution 
that may be driven by their comparatively large envelope 
densities. 



We also show in Figure [6] the c olors of two FHSC can- 
didates in_^erseus^_^er-Bolo 58 (jEnoch et al.ll20l'Ql ) and 
Bl-bS (jPezzuto et al.l[2012i) . In this diagram, the 70/24 
color of Per-Bolo 58 appears gen e rally consistent with 
that of a VeLLO, as lEnoch et all (|2010ii ) point out. As 
such, this source may be an extremely low-mass proto- 
star. On the other hand, the 70/24 color of Bl-bS is com- 
parable to the very reddest sources we find in Orion while 
the 70 /im flux is consistent with VeLLOs and fainter 
than the reddest sources in Orion by more than one order 
of magnitude. The faint but robust^etection of a 70 /im 
point-source by iPezzuto et al.l (j2012f ) may indeed point 
to the possible Class or VeLLO nature of Bl-bS. We 
do however note that Pezzuto et al. also detect a source 
with no 70 jam counterpart, Bl-bN, which may therefore 
represent a more robust FHSC candidate. Regardless of 
the elusive nature of FHSC candidates, when compar- 
ing ^oui^jiew_cajidlda^e^^ colors to FHSC models 
by ICommercon et aP (j2012l ) , we find that our sources do 
not appear to be consistent with predicted or expected 
FHSC colors, with the caveat that distinguishing FHSCs 
from VeLLOs with continuum observations alone is likely 
difficult. 

5. PACS BRIGHT RED SOURCES 

Up to this point, we have discussed two distinct and 
well-defined samples of sources in Orion: i) the sam- 
ple of candidate protostars identified with Herschel that 
have PACS 70 /im and 160 /im detections but MIPS 
24 jam magnitudes greater than 7.0 mag, ii) the sample 
of protostars that were reliably i dentified with Spitzer 
(24 /im mag:nitudes < 7.0 mag iMegeath et al.l 120121 : 
iKrvukova et al.|[2012f ). The protostar catalog target list 
used for the HOPS program consists mostly of the Spitzer 
identified protostars, but also contains some of the pre- 
viously known protostars with m(24) > 7.0 mag (Fischer 
et al., in preparation). 

In what follows, we focus our analysis on the Orion 
protostars that have 70/24 > 1.65. Of the 18 known 
protostars that satisfy this limit, 11 are identified with 
Herschel] hence this color regime is dominated by our 
newly identified protostars. Accordingly, Herschel has 
provided us for the first time with a far more complete 
sample of these red sources within the field of the HOPS 
survey. Given their red colors and their brightness in 
the PACS wavelength bands, we refer to this sample of 
protostars as PACS Bright Red sources, or PBRs. The 
coordinates, Spitzer photometry , and Herschel photom- 
etry of the sample are listed in Table HI while the APEX 
350 /im and 870 /im photometry are presented in Table [5l 
Since the APEX photometry are non-trivial to extract 
due in large part to contamination by cold surrounding 
material and, at 870 /im specifically, by the large beam 
size, we present three measures of the source flux, as de- 
scribed above. 

5.1. Observed properties PBRs 

As discussed above, we select 18 PBRs in Orion with 
observed 70/24 colors greater than 1.65. We show 4.5 /im 
to 870 /im images of 5 example PBRs in Figures [7] to [8] 
(see Appendix for the full sample images). Furthermore, 
in Figure [9] we show the full set of 18 PBRs observed 
SEDs from 24 /im to 870 jam. Inspection of the ob- 
served SEDs confirms that the PBRs sample is com- 
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Fig. 6. — Top: 70 /xm flux vs. 70 /xm to 24 /xm flux ratio for HOPS-detected Orion protostars. Open squares and open circles indicate the 
new candidate protostars with and without MIPS 24 /im detections, respectively. The color of the symbols indicates the flag values shown 
in Tableland discussed in § 3.2: here, green indicates sources that are reliable protostar candidates (flag = 1), orange indicates less reliable 
sources (flag = 2), red indicates extragalactic contamination (flag = 3), purple indicates other sources (flag = 4), and black indicates the 
single source without IRAC coverage. The solid vertical line indicates a 70 yum to 24 /xm flux ratio of 1.65, our PBRs selection criterion 
(see § 5). Triangles and x symbols indicate previously detected and characterized Spitzer protostars and faint protostellar candidates from 
[Megeath et al. (2012), respectively. The arrow indicates the extinction vector for a value of Ax = 6. For comparison, we also indicate 
as large diamonds the scaled observed measurements for some well-known Class sources, VeLLOs, and two FHSC candidates (see text). 
Bottom: 160 /im to 70 /im flux ratio vs. 70 /im to 24 /im flux ratio for HOPS-detected Orion protostars. The symbol and color type is 
the same as in the top panel. The reddest sources are distinguished from the bulk of Spitzer-identifLed HOPS sources (triangles) both by 
brighter 160/im fluxes and redder 70 /im to 24 /xm flux ratios. 
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Fig. 7. — 4^ x 4^ images of three PBR sources, 091015 and 091016 (top), and HOPS358 (bottom), at the indicated wavelengths, shown 
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posed of cold, envelope dominated sources with peak 
emission always located at A > 70 /im. In addition, 
the peak of the SEDs, and thus the temperatures, are 
well-constrained for all PBRs because we have obtained 
APEX sub-millimeter coverage for all sources. 

In Table [6l we present some basic properties of the 
PBRs. In particular, we find that 12/18 sources exhibit 
Spitzer 4.5 /im emission indicative of outflow activity. 
We also include some references to previous detections 
(see Appendix). Furthermore, 4/18 sources have signifi- 
cant levels of 4.5 /im emission that are indeed consistent 
with a high inclination. The majority of sources, how- 
ever, do not give clear indications of their inclinations 
at any observed wavelength, and therefore we cannot 
make any statements about their possible orientations 
based on their appearance in the images. We find in- 
dications from the 4.5 /im image morphology that two 
sources (HOPS341 and HOPS354) are binaries, while 
seven sources have a nearby source within 30". Two 
reside in more crowded regions, and seven sources ap- 
pear truly isolated. We find that a significant fraction 
(13/18) of sources appear to reside in filamentary regions, 
i.e., the extended 870 /im emission appears significantly 
elongated. 

The four sources with significant indications of a high 
inclination orientation are HOPS169, 302002, HOPS341, 
and HOPS354 (see Appendix for Figures [IH EH [23l 
and [25]). Inspection of their 4.5 /im images reveals that 
their outflows appear well collimated and relatively nar- 
row. Indeed, we might expect that sources that that 
have denser envelopes, and are therefore presumably 
younge r, may have more na rrow cavity opening angles 
(e.g., lArce fc SargentI l2006f ). As an additional check 
on our density analysis (see above), we use this inclina- 
tion information for an independent check of the enve- 
lope densities of these four sources. Despite the relatively 
sparsely sampled SEDs, we fix the inclination to 87° and 
fit the source SEDs. We find that even when we fix the 
model inclination to ^ = 87°, we still obtain envelope 
densities significantly above the pi value found in the 
previous section. 

5.2. Observational evolutionary diagnostics 

We calcula te Lhoi, Thoi (|Mvers fc Laddl I1993D . and 
Lsmm/Lboi (Andre et all I1993L 120001). The errors in Lboi 
and Tboi are derived with the same Monte Carlo method 
as described in § 6.2 for the modified black-body pa- 
rameters. We exclude the IRAC upper limits from this 
analysis; including these limits has an effect on our Lboi 
and Tboi estimates that is smaller than our estimated 
errors. We do, however, include the 24 /im upper limits; 
therefore the Lboi and Tboi values should be considered 
upper limits for sources not detected at this wavelength. 
Furthermore, we investigate the effect of applying an av- 
erage foreground reddening correction to all the new Her- 
schel candidate protostars. We find that dereddening the 
observed fluxes with extinction levels of Ay = 40 mag- 
nitudes has no effect on the derived parameters because 
the observed SEDs are extremely red and cold. 

5.3. Spatial distribution of the three samples 

We show the locations of the Herschel protostar candi- 
dates compared to the locations of the HOPS sample in 



Figure [TOl these positions are overlaid on the extinction 
map of Orion. It is immediately apparent that the spatial 
distribution of the new candidate protostars and PBRs 
is non-uniform. To investigate this distribution further, 
we show the relative fraction of new sources as a function 
of individual region in Table [71 The over-all number of 
new candidate protostars and PBRs is dominated by the 
Orion A cloud, and in particular L1641. This is not sur- 
prising since the LI 641 region is quite large and contains 
more protostars compared to other Orion regions. The 
fractions of new candidate protostars and PBRs com- 
pared to the total number of HOPS and new candidate 
protostars is, however, 2 times larger in Orion B. This re- 
sult is even more pronounced when we consider only the 
fractions of PBRs, with fractions that are more than 10 
times larger in Orion B. The NGC2068 (also containing 
the NGC2071 nebula) and NGC2024 (also containing the 
Horsehead or NGC2023 nebula) fields in Orion B have 
not only the largest fraction of new candidate protostars, 
but also of PBRs. While these numbers and fractions are 
subject to counting statistics and other possibly large 
sources of errors, the differences between Orion A and 
Orion B appear to be significant. 

About 5% of the combined protostars and candidate 
protostars in Orion are PBRs. If we consider the PBRs 
as representing a distinct phase in the evolution of a pro- 
tostar, and we assume a constant rate of star formation, 
the fraction the variation suggests that the protostars 
spend 5% of their lifetime in the PBRs phase (approx- 
imately 25,000 years with the 0.5 Myr protostellar life- 
time of lEvans eTal] [20091 ). averaging over all Orion re- 
gions. However, the assumed the duration of the PBRs 
phase would vary greatly with location, from 5, 000 years 
in the Orion A cloud to 80,000 years in the Orion B 
cloud. There are two alternative explanations. First, 
there might be environmental reasons which would favor 
the formations of PBRs, or perhaps extend the duration 
of the PBRs phase, in the Orion B cloud. Second, the 
ages of all the protostars in the Orion B cloud may be sys- 
tematically younger than those in the Orion A cloud. In 
this case, the regions containing the PBRs in the Orion B 
could be undergoing very recent bursts of star formation. 

Studies of pre-main sequence stars in the Orion 
molecular clouds show little evidence for significant 
ag:e differences bet ween the Orion A and B clouds. 
"Flahertv fc M uzerollei ((2008) determined an a ge of ^ 
2 Myr for NGC206 8 and NGC2071, while ReggianTeO] 
(j2011[ ). [Hsu et al.[ (12012) and lDa Rio et al.[ |2012) deter- 
mine ages for the ONC and LI 641 of ^2 — 3 Myr. How- 
ever, most of the protostars associated the NGC2068 and 
NGC2071 regions are outside the clusters of pre-main 
sequence stars a nd in dense filaments gas neighboring 
these clusters (Motte et al.ll2001ii) . A number of PBRs are 
in the LBS 23 clump (directly south of NGC2068) and 
in the NGC2023 clump (in the NGC2024 field); these 
are tw o of the 5 most rn assive, dense clumps found in 
Orion (iLada et al.[[l991al ). Compared to the other mas- 
sive clumps, both of these regions have 1/10 the num- 
bers of young st ars per unit gas mass and he nce may 
be quite young (jLada et al.[[1991b[ : [Ladal[1992[ ). Alter- 
natively, the gas in the LBS 23 and NGC2023 clumps 
may have dense gas filling fact ors that are much higher 
than the other massive clumps (jLada et al.[[l997l ): hence, 
they sources in these regions may be forming in a very 
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Fig. 9. — SEDs of the 18 PBRs are shown. The errors are smaller than the symbol size except for the 350 and 870 fim. data points. The 
black curve shows the modified black-body fit to the observed SED, with the indicated best fit parameters (see § 6.2). Note that PBRs 
19003 (shown in the top right) is located in a complex field where the photometry may be strongly contaminated. 
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different birth environment. The regions bordering the 
southern rim of the NGC2068 nebula and the northern 
rim of the NGC2Q71 nebula are also rich in protostars 
(see iMegeath et all 12012). Thus, the PBRs are found 
concentrated in sub-regions which may indeed be quite 
young. We will investigate these possibilities in future 
work. 

6. DETERMINING THE PHYSICAL PROPERTIES OF PBRS 
THROUGH MODELS 

In this section we describe our analysis of the physi- 
cal properties of the PBRs as inferred from their colors 
and SEDs. We first compare the 70/24 colors of the 
PBRs sample with those derived from a grid of mod- 
els which adopt the solution f or a rotating envelop e un- 
dergoing collapse fUh ichlllQTGUTerebev et al.lll98^ with 
outflow caviti es along the rotation axis of the envelope 
([Whitney et al. 2003a). This analysis puts a constraint 
on the minimum inner density of the protostellar enve- 
lope. Next, we compare the obser ved SEDs to m odel 
SEDs generated using the Hyperion (l Robitaill^l2011 [) ra- 
diative transfer code. This set of models assumes radial 
power-law gradients consistent with either a collapsing 
core with a constant infall rate and a static isothermal 
core. The models also encompass various combinations 
of internal and external heating. Given the prohibitively 
large computational time needed to explore the full range 
of parameters space using radiative transfer models, and 
given our inability to distinguish between models purely 
from five to six photometry points, we do not provide 
individual model fits to each protostars. Instead we fit a 
single temperature modified black-body function to the 
observed SEDs at 70 /im and longer wavelengths. The 
modified black-body fits provide luminosities and an ini- 
tial characterization of the envelope masses of the PBRs 
sample. 

6.1. Axisymmetric models: interpreting the 70/24 color 

We begin our ana lysis by using a simplified version of 
the lAli et all ()2010l ) protostellar envelope model grid to 
predict observed fiuxes and colors for comparison with 
our PBRs. The density distribution of these models is 
that of t he co l lapse of a spherical cloud in uniform ro- 
tation ([Ul rich"1976'), which is the inner region of the 
iTerebev e t al. (1984) model of the collapse of the slowly- 
rotating isothermal sphere. This model is then modi- 
fied by the inclusion of outfiow cavities of various shapes 
(jWhitnev et al.l [2003al lbl). This schematic model enve- 
lope captures the dependence of the short wavelength 
(24 /im and 70 /im) fiuxes on inclination due to rotation 
and bipolar cavities. 

The model fiuxes depend upon the mass infall rate, the 
angular momentum of the mass currently falling in, the 
outfiow cavity properties, the inclination of the rotation 
axis to the line of sight, and the luminosity of the cen- 
tral source, as well as the assumed dust properties. The 
thermal emission of the dusty envelope does not depend 
directly on the mass infall rate but instead on the density 
of the envelope. The model assumes free-fall at a con- 
stant rate, which results in a density profile with shape 
p oc r~^/^ (jTerebev et al.l 1198^ . The overah scaling of 
the density is characterized by pi, the density at 1 AU 



in the limit of no rotation: 

= l^ lO-eMeyr-i j^aSA^j ' 

The envelope mass infall Menv rate is related to p\ via 
the free-fall velocity, which in turn depends upon the 
unknown central mass M*. The actual model density 
structure departs from r~^/^ on small scales because of 
the angular momentum of the infalling material. This 
enters into the model through the parameter i^disk, the 
outer disk radius at which infalling material currently 
lands (see Ulrich 1976). 

The rotation leads to a significant dependence of the 
SED on the inc lination of the rotation axis relative to 
the line of sight (|Kenvon et al.| [r993) . This dependence is 
significantly enhanced by the inclusion of outfiow cavities 
(Whitney et al. 2003a) which are assumed to be aligned 
along the rotation axis. Finally, the overall shape of the 
SED is only weakl y affected by the lu minosity of the 
central source (L*; iKenvon et al.l [T993[ ) and so this is 
easily scaled. 

To roughly compare observed PBRs colors with those 
predicted by our model grid in Figure [11] we show the ef- 
fects of varying the model inclination, envelope density, 
and cavity opening angle on the 70/24 color and 70 /im 
fiux. As stated above, t hese mode l tracks are based on a 
simplified version of the lAli et al.l ([2010) model grid; we 
refer the reader to that publication for details. In brief, 
the model tracks that we consider here have the same 
fixed parameters as those listed in Table 1 of lAli et al] 
(2010), including a fixed central mass of 0.5 M©; in 
addition we have fixed the cavity shape exponent to 
a value of h = 1.5, and the envelope outer radius to 
^env,max = 1 X 10^ AU. We havc, however, expanded the 
envelope infall rate grid relative to the lAli et al.l (|2010l) 
grid to larger values (up to Menv = 10~^ yr~^ on a 
pseudo-logarithmic grid) and included a model with no 
envelope. As described above, for our model grid we as- 
sume that the envelope density falls off as p(r) ex r~'^l^ . 
In addition, while our model grid is parametrized in 
terms of Menv, with a fixed central mass of =0.5 M©, 
we will refer to p\ throughout this section. 

We note that the assumed central masses of protostel- 
lar sources re main largely unconstrained observationally 
(however see iTobin et al.| [2012). but may be lower than 
our assumed value. The effect of a lower central mass 
would be to lower Menv for a given value of p\. For 
example, if the central mass is 0.2 M©, the infall rates 
reported for our models would be reduced by a factor 
of 0.6. The assumed central mass, however, does not 
change the value of p\ corresponding to a given model 
SED. Furthermore, our small assumed disk radius does 
not strongly affect the trends shown here. 

With this model grid we isolate the effects, albeit in 
a simplified fashion, of varying the model inclination 
(viewing angle to the protostar), envelope density (pi), 
and cavity opening angle (Qq) on the 70 /im fiuxes and 
70/24 colors. In Figure fTH we show model tracks through 
70 /im fiux vs. 70/24 color space for high inclination (87°) 
viewing angles as a function of both and envelope den- 
sity (pi). By analyzing only the high inclination models 
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Fig. 10. — The distribution of red protostars overlaid on the extinction map of Orion. The dark grey contour shows the HOPS PACS 
coverage and the white contour shows the LADOGA 870 fim coverage. The triangles (light grey) indicate the positions of the HOPS 
protostars while the black points (squares and circles) indicate the positions of the new candidate protostars. The PBRs are highlighted in 
fuchsia and the properties of four selected PBRs are shown along with their corresponding 1^ x 1^ HOPS 70 fim images. 
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we obtain a lower limit on the envelope density required 
to reproduce the red 70/24 colors. 

We assume a fixed disk accretion rate of Mdisk = 1.0 x 
10~^ M0 yr~-^, a disk outer radius of Rdisk.max = 5 AU, 
and a central source luminosity of 1 Lq. Increasing 

Mdisk and Rdisk,max drives the models to both bluer 
70/24 colors and brighter 70 /im fluxes, while decreasing 
the inclination drives the models to bluer 70/24 colors. 
Therefore, we find that no model in our grid can ex- 
plain 70/24 colors with envelope densities less than log 
y9i/(gcm~^) ~ —13.4 (or equivalently an envelope infall 
rate of Menv = 5.4 x 10~^ M© yr"^). While models 
with larger values of pi and other combinations of pa- 
rameters can be found for bluer 70/24 colors, this anal- 
ysis sets an approximate lower limit on the expected 
envelope densities of sources with 70/24 > 1.65 of log 
Pi/(gcm~^) > —13.4. That is, high source inclinations 
alone cannot explain the red 70/24 colors of the PBRs, 
which also require dense envelopes. 

We expect that this value of pi sh o uld li e well within 
the Class range; e.g., iFurlan et al.l ()2008l ) found for a 
sample of 22 Class I sources in Taurus log pi/(gcm~^) < 
— 13.2, while over half had pi values that are lower than 
los Pi/(gcm~^) — —13.4. Note that the Fur Ian et al. 
([2008) data-set included far better sampling of the 
source SEDs and, most critically, Spitzer IRS spec- 
troscopy, allowing for a more robust estimate of the pa- 
rameters of sources in their sample. In contrast, our 
SEDs generally are envelope dominated with few sources 
having robust detections shortward of 24 pm. Therefore, 
any estimate of the value of the envelope density will 
be necessarily imprecise and have a back of the envelope 
character. That said, our comparison with the model 
grid and previous derived values of the envelope density 
leads us to conclude that the 70/24 > 1.65 color cut, 
while not uniformly selecting a unique envelope density 
threshold, will preferentially select Class sources with 
log pi/(gcm~^) > —13.4, irrespective of source inclina- 
tion. 

Having demonstrated that the very red 70/24 colors 
require large values of pi, as expected for SEDs peaking 
at long wavelengths, we proceed by analyzing the PBRs 
in the context of more simplified models. In what fol- 
lows, we carry out two independent analyses of the source 
SEDs: a qualitative model image comparison using_JT) 
models generated with the Hyperion (Robitaille 2 011f ) 
radiative transfer code, and modified black-body fitting 
to the observed SEDs at 70 pm and longer wavelengths. 
We attempt to maintain consistency by using similar en- 
velope dust models throughout the analysis presented 
here. For the Hyp erion model image analysis we use the 
lOrmel et al.l ([2011) opacities. These opacities are similar 
to the commonly assumed Ossenkopf & Henning (1994) 
("0H5"; see below) opacities but include both the scat- 
tering and absorption components at short wavelengths, 
needed for radiative transfer calculations. Specifically, 
we use the "icsgra2" lOrm el et al. (2011) model opaci- 
ties, which include icy silicates and bare graphites, with 
a coagulation ti me of 0.1 Myr. We as sume the dust 
model from Ossenkopf fc HenningI ()1994f ) for the long- 
wavelength modified black-body fits, specifically, their 
"0H5" opacities, corresponding to column 5 of their Ta- 
ble 1. These opacities reflect grains having thin ice man- 



tles with 10^ years of coagulation time at an assumed gas 
density of 10^ cm~^. 

6.2. Spherically symmetric models: Comparison to 
photometry derived from model images 

We use the Hyperion (Robitaille 2011) Monte Carlo ra- 
diative transfer code to investigate some basic properties 
of the PBRs. We run a series of spherically symmetric 
models under a range of very simple assumptions, and 
produce simulated images. As stated ab o ve, th e dust 
model we assume is that of lOrmel et al.l (poTT) (" ics- 
gra2"). To explore which model assumptions might be 
reasonable for analyzing the properties of the reddest 
sources in Orion, we investigate a few limiting model 
scenarios. We therefore disregard the details of the indi- 
vidual source SEDs of the entire red sample and qualita- 
tively compare the SEDs of the two sources with the 
smallest and largest values of Lboi and Tboi (sources 
091016 and HOPS358, respectively; see Table Ej) with 
photometry derived from model images. 

We generate a series of model images that fall into 
four classes: i) a starless core at a constant tempera- 
ture of 10 K (referred to here as the "core" model); ii) 
a starless core with an isotropic external radiation field 
("core+external"); iii) a core with an internal source 
("star" model); and finally, iv), a core with both an in- 
ternal source and an isotropic external radiation field 
("star+external"). For each of these classes of mod- 
els, we test a range of density profile shapes and den- 
sity normalizations. For the density profile shape we 
assume two values, a = 1.5,2.0, where a is the ra- 
dial density profile power law index: p{r) ex r~". For 
the absolute value of the (gas) density normalization at 
1 AU we assume 5 values: log (piAu/(^ cm~^)) = —10 
to —14, in steps of (5 log (piAu/(^ cm~^)) = 1.0. For 
models with external heating, the bolometric strength 
of the int erstellar radiation f ield (ISRF) is set to the 
value from iMat his et al.] (j 19831 ) at the solar neighborhood 
(47r = 0.0217 ergs/cm^/s). The spectrum of the radia- 
tion field is assumed to be that at the solar neighborhood 
from Porter fc Strong (2005), but reddened by Ay 10 
using the Ikim et al.l fl994) extinction law. The ISRF 
model includes contributions from the stellar, PAH, and 
FIR thermal emission. The inner radius of the core is set 
to the radius at which the dust sublimates, assuming a 
sublimation temperature of 1,600 K, while the outer ra- 
dius is set to Ipc. The central source is taken to have 
10 L0 and a spectrum given by a Planck function at the 
effective temperature of the Sun (5778 K). The choice 
of the stellar temperature is arbitrary, and is unimpor- 
tant for the modeling presented here, since all sources are 
deeply embedded and all stellar radiation is reprocessed 
— on ly the total bolometri c luminosity is important (see 
e.g., iJohnston et al.ll2012L for a discussion of the Rstar 
and Tstar degeneracy). 

The high levels of spatial filtering caused by our 
adopted aperture photometry scheme require us to as- 
sume such a central source luminosity to roughly match 
the flux levels in the observed SEDs. Furthermore, while 
a 1 pc sized envelope is larger than usually assumed, the 
high levels of spatial filtering inherent in the aperture 
photometry cause us to be insensitive to structure on 
scales larger than the assumed aperture sizes. 
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Fig. 11. — Left: Dependence of the model 70 /i,m to 24 fim flux ratio for high inclination orientations (87°) on envelope density (pi) 
for the 5 cavity opening angels in our grid. These model tracks assume an internal luminosity of 1.05 Lq. The 70/24 color selection is 
indicated with the vertical black line, while the corresponding log pi/(g cm~^) ~ —13.4 (corresponding to an envelope infall rate of Menv 
= 7.5 X 10~^) threshold is indicated with a horizontal black line. Right: Model 70 pm flux vs. 70 pm to 24 pm flux ratio for high inclination 
orientations (87°). For the model with a cavity opening angle of 6c = 25^ (H — symbols) and log pi / (g cm~^) ~ —13.4, we show the effect 
of of increasing the internal luminosity by two orders of magnitude with the dashed grey curve while the effect of decreasing the inclination 
to 18° is illustrated with the dot-dashed curve. All models with higher assumed values of M^isk and lower inclination will have envelopes 
that are denser than log pi/{gcm~^) ~ 13.3, the median value of pi in our model grid; sources with redder 70 to 24 pm. colors cannot be 
explained by envelopes that are less dense than this threshold. The light grey thick line corresponds to a 24 pm limit of 7 magnitudes, the 
imposed Spitzer protostar magnitude limit for identiflcation of protostars. 



The model images have a resolution of l"pix~^^ or 
420 AU at our assumed distance. We convolve these 
images with the azi muthally averaged PSFs provided by 
I Aniano et aH (j2Qll[ ). except in the case of the SABOCA 
350 fim and LABOCA 870 /im images. These wave- 
lengths are convolved with Gaussian PSFs with FWHMs 
equal to 7.4'' and 19'' respectively, i.e., the nominal beam 
sizes for our observations. All model image photometry is 
then performed on the convolved model images using the 
same aperture and sky annulus parameters as those ap- 
plied to the data. The use of such photometric aperture 
parameters can cause large amounts of spatial filtering 
due the small sizes of the apertures relative to the beam 
sizes and the extent of the core emission (see below). 

We show our extracted model SEDs in Figure [12] and 
a subset of the corresponding model images in Figure [T3l 
In general, we find that models without internal sources 
are very unlikely to match the observed PBRs proper- 
ties, on the basis that their SEDs are less luminous and 
peak at longer wavelengths than the observed SEDs for 
all density profile shapes that we assume; see Figure [121 
top panels. While the " core+external" models (top right 
panel) suffer from severe spatial filtering, it is unlikely 
that such models will well represent the data as these 
SEDs tend to also peak at longer wavelengths. For the 
two classes of models with internal sources we find bet- 
ter agreement with the data; see lower panels of Fig- 
ure [121 While steeper envelope profiles may imply some- 
what higher envelope densities, the range in plausible 
densities for a = 1.5 is logpi ^ —12 to —13, while for 
a = 2.0 values of log pi ~ —12 roughly agree with the 
shapes the observed SEDs. We note that we find very lit- 
tle difference between the "star+external" and the " star" 
model in the A < 160 /im regime, possibly indicating that 
our assumed ISRF strength relative to the assumed inter- 
nal source luminosity may be underestimated compared 



to what we might expect to find regions like Orion. 

6.3. Modified black-body fits to the PBRs 

An individual detailed fit to each protostar is be- 
yond the scope of this work considering the vast amount 
of parameter space needed to model protostellar SEDs 
(i.e., source luminosity, envelope density, envelope ro- 
tation, outfiow cavity geometry, external heating, outer 
envelope structure). Furthermore, unlike most of the 
HOPS protostars, whose properties can be constrained 
from a combination of far-IR photometry, 5 /im to 
40 / im Spitzer /IRS spectra, and Hubble near-IR imag- 
ing f Fische r et aTll2010L "2012\ the properties of the SEDs 
must currently be derived from 5 to 6 photometry points 
at long wavelength. While the above modeling and anal- 
ysis shows that the internal source is important, the 
longer wavelength fiuxes are probing the bulk of the 
envelope mass, expected to mostly be at a single tem- 
perature. Furthermore, for density profiles in the range 
of p{r) (X r~^/^ or r~^, as assumed above, we expect 
that most of the envelope will be located at large radii. 
We therefore perform modified black-body fits to the 
A > 70 /im SEDs hsted in Tables [4] and [5l For this anal- 
ysis we use the beam fiux measurements for the sub- 
millimeter 350 /im and 870 /im portion of the SED. The 
results of the analysis are presented in Table [8] and the 
model SEDs are plotted with the data in Figure [9l 

Before fitting the long-wavelength SEDs of the sources, 
we apply color corrections to the Herschel 70 /im, 
100 / zm, and 160 /im fiuxes. Following (^ Launhardt et al.1 
l2012f ) , these photometric color corrections have been de- 
rived iteratively from the slopes of the PACS SEDs, us- 
ing polynomial fits to the values in Table 2 of the PACS 
calibration release note "PACS Photometer Passbands 
and Colour Correction Factors for Various Source SEDs" 
from April 12, 2011. The color corrections for the APEX 
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Fig. 12. All panels show the SEDs of HOPS358 (blue; Lboi = 30.6 L© and Tboi = 44.2 K) and 091016 (red; Lboi = 0.65 and Tboi = 
29 K), the extrema of the Lboi and Tboi distributions for the reddest sources, compared to spherical models. The two line styles indicate the 
different assumptions for the density profile shape while the shade of grey indicates the different model density pi values. The four panels 
correspond to the four different models that we test (see text): Top left: the "core" model, a starless core; Top right: the "core+external" 
model, a starless core with external heating; Bottom left: the "star" model, a core with an internal source; and finally. Bottom right: the 
"star+external" model, consisting of a core with an internal source which is irradiated by an external radiation field. Severe spatial filtering 
due to our aperture photometry on the model images can be seen most notably in the two externally heated models (left column). 



data are assumed to be negligible. 
The form of the modified black-body function is given 

by 

5, = nB,(i.,Td)(l-e-^('')), (2) 

where ft is the solid angle of the emitting element, B^yiT^^) 
is the Planck function at a dust temperature Td, and t{v) 
is the optical depth at frequency v. Here, the optical 
depth is given by T(h>) = Nn mn RgJ^ ^(^), where A^h = 
2 X A^(H2) + A/'(E[) is the total hydrogen column density, 
mn in the p roton mass, Ku is the a ssume d dust opac- 
ity law from lOssenkopf fc HenningI (Il994f), and Rgd is 
the gas-to-dust ratio, assumed to be 110 (|Sodroski et ahl 
IT997I ). The best-fit total masses Mtot reported in Table [8] 
have been multiplied by an additional factor of 1.36 to ac- 
count for helium and metals. Furthermore, in Table [8] we 
also report the peak wavelength of the best-fit modified 
black-body model. Finally, we estimate Lgmm from the 
model SED, where Lgmm is integrated over A > 350 /im. 

If a given source SED has coverage over fewer than 
4 long-wavelength points, we do not fit a model to the 
SED. While all the PBRs sources satisfy this criterion. 



all the new candidate protostars and HOPS sources do 
not (see § 5.2) and are therefore not fitted. The errors on 
Td, Mtot, and the thermal component of the luminosity 
(Lmbb) are estimated through a straight-forward Monte 
Carlo methocQ- For each source we generate 2000 syn- 
thetic SEDs drawn from a normal density with mean and 
standard deviation equal to those of the measured SED 
at each wavelength. We then fit each synthetic SED. The 
reported error is equal to the standard deviation of the 
resulting distribution of each parameter. These errors 
do not include systematics introduced by, e.g., our dust 
model assumption or variation in the gas-to-dust ratio. 

We show the modified black-body fit results, along 
with the SEDs of the PBRs, in Figure M The resulting 
best-fit mass, luminosity, and temperature is also indi- 
cated for each source. The model fits the data surprising 
well considering that significant temperature gradients in 
the envelope are expected. Furthermore, in all cases the 

"Offered the choice between mastery of a five-foot shelf of 
analytical statistics books and middling ability at performing sta- 
tistical Monte Carlo simulations, we would surely choose to have 
the latter skill." Press, 1993, Numerical Recipes, page 686. 
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Fig. 13. — Convolved model images at the indicated wavelengths, shown on a log scale with the same minimum and maximum flux levels 
for all panels to illustrate the shape of the model SEDs. Each image is 200^^ (or 84000 AU) on a side. The top row shows the "core" model 
images, with log(pi) = —12 and a density proflle shape of p oc r"^ *-*, corresponding to the top left panel in Figure [12] The bottom row 
shows the "star+external" model, with the same density paramaters as the top row, corresponding to the bottom right panel in Figure fT2] 



24 jam point, when detected, has a much higher flux level 
than the modified black-body model. We interpret this 
discrepancy as strong evidence for internal heating by a 
protostar. 

Excluding the 70 /im point and fitting only the A > 
100 /im SED has a minor effect on the resulting param- 
eter values. Without 70 //m, the masses systematically 
increase by 40%, the temperatures decrease by 5%, and 
the luminosities decrease by 7%. This small effect may 
be understood by the fact that the 100 /im fluxes are 
well-correlated with the 70 /im fluxes for this sample, 
tracing similar material near the protostars. The tem- 
peratures at 100 /im and 70 /im are not dramatically 
different, and most likely both points are dominated by 
optical-depth effects such that the r = 2/3 surface is 
not s ignificantly differe nt between the two wavelengths. 
From iHartmannl (p009), the radius of the r = 2/3 surface 
can be roughly approximated as rx oc a^^; this relation 
implies that rro/rioo ~ 2.5. 

Our best-fit modified black-body model always under- 
estimates the observed 870 /im flux of all sources. The 
model sub-millimeter SEDs are always bluer than the 
observed SEDs. We find that the discrepancy is at the 
0.83 ±0.26 Jy level (or a factor of ~ 3) excess), where the 
error bar represents the standard deviation in the resid- 
ual distribution. It is likely that this discrepancy is dom- 
inated by the larger beam size of the 870 /im observations 
which has the effect of mixing the source flux with that of 
the surrounding cold and possibly high-column environ- 
ment. Contamination to the 870 jam fl ux by disk emis- 
sion m ay also increase this discrepancy. Ij0rgensen et al.l 
(pOOO) find average disk masses of ^ 0.13 (with a 
large scatter) across their sample of Class sources. The 
sources in their sample that are comparable to our PBRs, 
however, are those with the lowest values of T^^/. For ref- 
erence, they found that IRAS4A1, with a T^oi = 43 K, 
has the largest disk mass of 0.46 in their sample; on 
the other hand, LI 157, with a similar Ti^oi = 42 K, has a 
disk mass about a factor of 4 smaller. We estimate that 
a 30 K disk of 0.5 Mq would cont ribute ^ 0.6 Jy to the 
beam fiux at 870 /im (assuming lOssenkopf fc HenningI 



(|1994l ) dust opacities, as above). Therefore, disk emis- 
sion could indeed contribute to the observed 870 jam fiux 
but further detailed observations at high resolution are 
needed to disentangle the envelope component from the 
possible disk emission. Another possible source of am- 
biguity in interpreting the 870 /im fiux discrepancy is 
the model dust opacity assumption. Furthermore, we do 
not find an 870 /im discrepancy in the analysis of model 
images presented in the previous section. This indicates 
that large disk masses may not be necessary to explain 
the sub-millimeter fiuxes. We therefore emphasize that 
the disk masses inferred here from the 870 /im excess 
should be regarded only as upper limits; further detailed 
investigation into the disk properties of our sources is 
deferred to future work. 

Independent of these issues, it is clear that a more ac- 
curate treatment of the data would require all images to 
be convolved to a matched resolution; however, this ap- 
proach would have the effect of causing non-detections 
for a majority of sources at the shorter wavelengths due 
to the relatively large limiting beam size of our data-set 
(~ 19'' at 870 /im). Homogeneously extracted SEDs are 
therefore not feasible for this data-set as a whole. Nev- 
ertheless, we test the effects of convolving the the data 
before extracting the SEDs. We choose PBRs 119019 as 
a test source because it is isolated and has approximately 
median values for the best-fit modified black-body tem- 
perature, luminosity, and mass. Ignoring the 870 /im 
data, this source is clearly detected at 70 /im, 100 /im, 
160 /im, and 350 /im. For these four wavelengths, the 
largest beam size of ~ 12'' corresponds to the 160 /im 
data. We therefore convolve the 70 /im, 100 /im, and 
350 /im data to a resolution matching the 160 /im obser- 
vations and extract a beam-smoothed SED. We then fit 
this SED in the same way as described above. Compared 
to the non-convolved SED modified black-body fitting 
results, we find that the temperature decreases by ~ 4%, 
the luminosity increases by ~ 10%, and the the mass in 
the thermal component increases by ~ 30%. These sys- 
tematic shifts are similar to but somewhat larger than 
the errors quoted in Table [8] (~ 2% on the temperature. 
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~ 10% on the luminosity, and ~ 30% on the mass). 
We note, however, that the errors quoted in Table [8] are 
purely random and do not include any systematic compo- 
nent. We therefore conclude that extracting SEDs from 
images matched to a resolution of ~ 12" will not greatly 
affect our results. 

Modified black-body fits provide a somewhat limited 
means of analysis of our sources since the model assumes 
a single temperature and density along the line of sight 
for the emitting material. We expect that the assump- 
tion of a single line-of-sight temperature will caus e an 
unde r estimate of the source masses (e,g., Niel bock et al.l 
I2Q12I : iLaunhardt et al.|[2QT2l ). However, radiative trans- 
fer models have large ambiguities in the assumed source 
temperature and density structure, leading to mass esti- 
mates that strongly model-dependent. Furthermore, the 
dust law that is assumed will introduce significant uncer- 
tainties into the derived masses, irrespective of the analy- 
sis method that is implemented. For example, the masses 
listed in Table [8] increase by a factor of ^ on average 
when we assume Draine Led ([1984) Ry = 3.1 ISM- 
like dust. These issues indicate that the masses derived 
here represent lower limits to the true envelope masses. 
Nevertheless, we consider the modified black-body fits 
to the measured photometry to provide the most robust 
estimates of the mass that we currently have. 

We note that with only 5 SED fiux points at best, fit- 
ting a multiple component (modified) black-body model 
cannot be justified. Since most of the mass is located 
at relatively large scales and expected to have cold tem- 
peratures, excluding the warmer shorter wavelength data 
arising from inner material will not significantly increase 
the masses we derive. The modified black-body fits thus 
to provide an approximate measurement of the optical- 
depth averaged gross properties of the envelopes being 
investigated. These issues point to the need for a more 
sophisticated modeling approach that will be carried out 
in future work. 

7. DISCUSSION 

As seen in Figure [TTJ the observed 70/24 colors of a 
protostar can be driven towards redder values through 
various strongly degenerate parameters. For example, 
the total column of material along the line-of-sight 
(LOS) towards a given protostar can have multiple con- 
tributions: the attenuation of the mid-IR emission by 
dense foreground material, the density of the envelope, 
the amount of envelope fiattening, the opening angle of 
the outfiow cavity, and the source inclination. Further- 
more, the assumed model central protostar mass remains 
largely unconstrained by observations to date and can af- 
fect the interpretation of the 70/24 colors. 

Foreground extinction can have various contributions, 
such as intervening dust between the observer and the 
cloud and dense material associated with the cloud itself, 
such as filamentary material. Of these two components, 
the first is expected to be relatively small, while the latter 
can be expected to vary from source to source by rela- 
tively large amounts, with a corresponding effect on the 
observed colors. For example, some PBRs are located in 
filamentary regions (e.g.. Figures [7] and [8]) , while others 
appear more isolated (e.g., see Appendix for Figure [2Q|). 
We have estimated the effects of foreground extinction 
levels up to a level of Ay = 40 mag, and find that the 



values of Lboi, Lgmm, and Tboi not significantly affected. 

On the other hand, the effects of source inclination are 
not as straightforward to assess. When considering the 
presence of fiattened rotating envelopes, disks, and out- 
fiow cavities, the source inclination will have a large effect 
on the observed source SED, as illustrate d by the model 
tracks shown in Figure ITTI fsee also, e.g., IWhitnev et al.l 
120033 iDe Buizer et al.l[2005l : lOffner et al.l l2012'). 

Therefore, it appears that the very red PBRs can be 
explained by multiple effects that all result in increas- 
ingly red observed 70/24 colors. These very red colors 
may be driven by elevated envelope densities (or equiva- 
lent ly, Menv), high source inclinations, or elevated levels 
of extinction associated with structures larger than the 
envelope-protostar system. The current data and SED 
coverage do not allow us to break these degeneracies con- 
clusively. Furthermore, we consider it likely that the red 
observed colors are not driven any single cause, but in- 
stead are the result of several. 

We have, however, designed our PBRs selection to find 
the densest envelopes in Orion (c.f.. Figure [TT]). Further- 
more, the effect of external foreground extinction is not 
expected to be large at these long wavelengths, even with 
elevated levels of material along the LOS (see above). In- 
deed, we have also shown that the PBRs require a central 
heating source, indicating that the detection of a 70 /im 
point source drives the interpretation of the sample as 
Class sources, irrespective of source inclination. How- 
ever, we note that if the central masses are significantly 
different than the assumed value of 0.5 M©, then for a 
fixed reference envelope density the inferred envelope in- 
fall rates need to be scaled accordingly (see Equation 1). 

To further investigate the evolutionary state of the 
PBRs, in Figure [TH we show we the values of Lboi, 
Tboi, and Lgmm/Lboifor the entire sample of new can- 
didate protostars and the previously identified Spitzer 
HOPS sample (Fischer et al. in preparation). In the 
left panel, we show Lboi vs. Tboi for the entire sam- 
ple of new protostar candidates, including those fiagged 
as possible extra-galactic contamination. We also show 
the four reference Class sources presented in Table [H 
The PBRs sample in particular, and the entire sample of 
new candidate protostars, are generally clustered around 
low Tboi values. Ignoring inclination degeneracies and 
other considerations, these low Tboi values indicate that 
the PBRs sample is indeed composed of young Class 
sources. In the right panel, we show Tboi vs. Lgmm/Lboi 
for the sources with sufficient coverage to estimate Lgmm 
(see § 6.2). The PBRs, as expected if the sample can 
be explained as Class sources, cluster around larger 
valu es of Lsmm/Lhoi co mpared to the rest of the sam- 
ple. I Andre et al.l (|2000f ) proposed the Lgmm/Lboi > 0.5% 
threshold for Class sources, and all but one of the new 
candidate protostars for which we can estimate Lgmm fall 
into this category. Irrespective of the evolutionary in- 
dicator that is chosen (Tboi or Lgmm/Lboi ), all of the 
new candidate protostars in both the reliable and lower 
probability categories (green and yellow points, respec- 
tively), would be considered of Class status. Finally, 
while the PBRs 70/24 > 1.65 color criterion causes some 
sources with very low values of Tboi and very high values 
of Lgmm /Lboi to be missed, the color selection is able to 
capture the vast majority of the most extreme Class 
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Fig. 14. — Left: Bolometric luminosity vs. bolometric temperature for the new candidate protostars (H — symbols) and the HOPS protostars 
(light-grey triangles). The squares indicate the PBRs sample, drawn from both the new candidate protostar sample and the HOPS protostar 
sample. The dashed lines indicate the canonical Tboi divisions between protostellar Classes. Right: Lsmm/Lboi vs bolometric temperature 
for the subset of sources shown in the left panel for which we have sufficient submillimeter coverage to estimate Lsmm (see text). The 
horizontal dotted line indicates the Andre et al.. (,2000 ) proposed Lgmm/Lbol > 0.5% Class threshold. 
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sources at the extrema of the Lgmm/Lboi and Tboi distri- 
butions. 

This evidence strongly supports the interpretation of 
the PBRs (and indeed ah the sources classified as reli- 
able protostellar candidates) as very dense Class pro- 
tostars, irrespective of the source inclination. On the 
other hand, the new candidate protostar sample, taken 
as a whole, may be explained by a combination of the 
effects described above: high inclination, high densities, 
and extreme values of foreground extinction, along with 
elevated levels of extragalactic contamination. Of par- 
ticular interest is the possibility that some of the sources 
classified as low probability protostars at low Lboi values 
may be confirme d as bona fide protost ars with future 
observations (see lOffner fc McKeell2Qlll for a detailed 
discussion of the significance of such sources); this will be 
investigated in future work (Stutz et al. in preparation). 

To definitely measure the inclinations of our sample 
of sources and to therefore determine envelope densi- 
ties more accurately we require millimeter line emission 
maps at high resolution from e.g., ALMA, along with 
single dish observations of high density tracers. A vig- 
orous follow-up campaign is therefore underway to more 
firmly place these protostars within the context of star 
formation in the Orion clouds. We are observing am- 
monia spectra toward the full sample of protostar can- 
didates to verify the presence of dense molecular gas 
and determine kinetic temperatures. A Herschel PACS 
range spectroscopy program toward 8 PBRs will char- 
acterize the energetics of outflows and UV heating on 
small-scales. The outflows (CO), dust continuum, and 
surrounding dense molecular gas (N2H+) are being ob- 
served in the millimeter to determine the source inclina- 
tions, outflow opening angles, inner envelope properties, 
and the kinematics of the larger-scale dense gas. 

8. CONCLUSIONS 

We have discovered a sample of 55 new candidate pro- 
tostars in Orion with Herschel, as part of the HOPS 
Open Time Key Programme scan-map observations at 
70 /im and 160 /im. We conclude that: 

• The new candidate protostars are either very faint or 
undetected at Spitzer wavelengths. We find 34 sources 
with 24 /im magnitudes that are greater than 7.0 and 21 
sources that are undetected in the MIPS 24 /im band. 

• We analyze the IRAC colors and the broad shape of the 
SEDs between 3.6 jam and 160 /im. Based on this anal- 
ysis, we classify the sample as follows: 27% (15 sources) 
are considered reliable protostars, 18% (10 sources) are 
considered lower probability protostars, 47% (26 sources) 
are classified as extragalactic contamination, including 
AGN, 3 sources have IRAC colors consistent with stel- 
lar photospheres but Herschel and APEX SEDs consis- 
tent with cold dust emission, and 1 source does not have 
IRAC coverage. We find that the subset of sources with- 
out 24 /im detections dominates the number of sources 
categorized as most reliable protostellar candidates by 
a factor of ~ 3, suggesting that sources with no short- 
wavelength detections and only PACS 70 /im and longer 
wavelength detections are much less likely to be of extra- 
galactic origin. 

• We combine the new protostar candidate sample with 
the previously identified Spitzer HOPS sample and find 
that 18 sources have 70/24 colors greater than 1.65. 



These are the reddest protostars known in Orion, 11 of 
which are newly identified Herschel sources in the re- 
liable protostar category listed above. We name these 
sources "PACS Bright Red sources", or PBRs. Com- 
pared to the other protostars in the HOPS fields, the 
PBRs populate the extrema in the distributions of stan- 
dard evolutionary diagnostics, having both the largest 
Lsmm/Lboi ratios and lowest Tboi values. The PBR 
source SEDs and peak SED wavelengths are consistent 
with the hypothesis that the PBRs do indeed represent a 
population composed of Class sources with the densest 
envelopes in Orion. 

• A comparison to radiative transfer models of rotat- 
ing, collapsing protostellar envelopes with outflow cavi- 
ties show that the 70/24 > 1.65 color limit selects sources 
with envelope densities with log pi/(g cm~^) > —13.4, ir- 
respective of inclination effects. While the 70/24 > 1.65 
color selects sources with dense envelopes, this color cri- 
terion does not find all dense sources above the thresh- 
old of of pi/(gcm~^) = —13.4. Therefore this selection 
should be used in conjunction with other evolutionary in- 
dicators, namely millimeter and sub-millimeter measure- 
ments, to determine the nature of the observed sources. 

• Our modeling of the PBRs SEDs reveals that these 
sources are not consistent with being externally heated 
starless cores; the presence of a 70 /im point source re- 
quires that the sample be interpreted as dense envelopes 
containing embedded protostars. 

• The fraction of known protostars that are PBRs varies 
from 1% in the Orion A cloud to 17% in the Orion B 
cloud, with an an average fraction over the Orion com- 
plex of 5%. These numbers suggest that if the PBRs 
represent a distinct phase in protostellar evolutions, pro- 
tostars spend on average 5% of their lifetime in the PBRs 
phase. Most of the PBRs in the Orion B cloud are con- 
centrated in dense gas near the NGC2068, NGC2071, 
and NGC2023 nebulae. These regions of dense gas are 
also known for a lack of more evolved pre-main sequence 
stars. The high percentage of PBRs in Orion B suggests 
either that the regions containing the PBRs may cur- 
rently be undergoing more vigorous star-formation than 
other regions of Orion, or that the PBRs lifetime is longer 
in these regions. 

• The sources with 70/24 < 1.65 colors and faint 70 /im 
fluxes must be confirmed as either protostellar or con- 
tamination sources before their significance can be as- 
sessed. If confirmed as the former, however, this sample 
would constitute an important population of very low- 
luminosity, cold protostars previously unobserved, most 
interesting from the point of view of constraining the 
faint end of the luminosity function of protostars. 

• We expect that a comparable number of very red proto- 
stars will be found with Herschel in more near-by star- 
forming regions. We caution that a careful treatment 
of possible extragalactic contamination must be imple- 
mented to understand the broader significance of such 
sources. 
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APPENDIX 

PREVIOUS PBRS DETECTIONS 

We searched the SIMBAD and VizieR (|Ochsenbein et alJ l2QQQ) services for previous identifications of the PBRs. We 
restricted our search to a radius of 20'' from the 70 /xm source coordinates. The results are summarized in Table [6l 
This list is likely incomplete and is intended to provide a resource and rough guide to some of the previous detections 
of these sources. 

— 061012, 119019, and 097002: None found. 

— H OPS169: Known protostar and outflow bipolar outflow V380 Ori NE (e.g., iDavis et all 120001. "Stanke et_aD 
(2002) detected source 59 in their 2.12 /im catalog of Orion A, offset by from our source coor dinates. Davis et al. 
( 20091 ) also detected an outflow about 2" away from our coordinates. iNutter fc Ward- Thompson! (|2007[ ) classified this 
source as hosting a young stellar object (YSO) and measured an 850 /im-fiux-derived mass of 2.8 (assuming a 
temperature of 20 K and a distance of 400 pc). 

— 019003: Since this source is located in the very crowded and comp lex ONC filament, it is often not clear which 
previous identifications may be associated with it in particular. Tsujimoto et al.l (|2003.) listed a near-IR source located 
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3.35'' away from our coordinates. iNutter fc Ward-ThompsonI (|2QQ7[ ) detected a source offset by 8.4''; they classified it 
as hosting a YSO and measure a mass o f 6.2 M^:^. While it is not cle ar if the detection is associated with the PBRs we 
consider it probable. iChini etld] (|1997[ ) and iNielbock et ahl (|2QQ3[ ) likely detected this source in their 1.3 mm maps 
near FIRla in OMC-2, although their beam sizes were too large identify the source unambi guously. 

— 082005: This source was p reviously classified as a starless core (e.g., [Johnstone & Ballv 2006 
INutter & Ward-Thompson 2QQ1 fMookeriea et a l. 2009). The lNutter Ward-ThompsonI (l2007H ocation of the source 
is about 8" away from our source with an associated mass of 3.3 M0. iMookerjea et al.l ()2009l ) measured a mass of 
4.3 Mq for this source. 

— HOPS372 and 082012: These sources have coordinates that are significantly offset from the 
INutter & Ward-Thompson| ()2007f ) coordinates; the two so urces are unresolved and lie about 13" away from their 
catalog entree, with a mass of 12.7 M©. The IMookerjea et al. (2009) analysis derived a mass of 1 for 082012 
(MMl) and 7.4 M^;^ for HOPS372 (MM2), from modified black-body fits to long-wavelength SEDs. 

— 090003: Miettinen et al.l ()2009,) detected this source (SMM3) in their 850 /im map, and concluded that it is a 
promising Class cand idat e based on the shape o f the SED. They calculated a mass of about 7.5 for this source. 
iMiettinen et al.l ()2010l ) and iMiettinen et al.l (|2012f ) also observed this source in various molecular line transitions and 

with SABOCA at 350 /xm; their measured flux for this source is S^^q^ = 3.63 Jy beam~^. We include their flux 
measurement and the 350 um map in our analys is. 

— HOPS358: .Strom et al.l ([T976[ ).ISt rom et all ()1986f ) , and (i Reipurth et al.|[l999f ) detected the Her big-Haro complex 
HH24-26, with HH25 located about 10" away from HO PS358. This source is also included in the IWu et al.l l\2004! ) 
high velocity outflow catalog. INutter fc Ward-ThompsonI (j2007l ) classified a nearby source (11" away) as a starless 
core with a mass of 6.3 Mq. 

— 091015 and 091016: These sources were detected by iLis et al.1 (e.g., 1999) at 1.3 mm and 350 /im (as 
sources 5 and 6). The reported masses are ^ 2 Mq for each source. Both sources are classified as starless by 
INutter fc Ward-ThompsonI l\2007\ ), with catalog masses of 0.8 Mq (0 91015) and 1.3 M^ (091016 ). 

— HOPS373 and 093005: Both sources were classified as s tarless bv INutter fc Ward-ThompsonI ()2007[ ). with reported 
masses of 4.2 M© (HOPS373) and 3.9 M© (093005). Motte et al.l (j2001f ) observed this region at 450 and 850 /am and 
also classified 093005 as a starless core, based on the lac k of evidence for an embedded source. Source HOPS373, 
on the other hand, is known to be driv ing a CO outflow (Gibb fc Lit tlell2000f). and therefore has been classified as a 
candidate Cl ass source. Furthermore. iHaschick et al.l ()1983f ) identified a water maser near the HOPS373 location. 

— 302002: INutter fc Ward-ThompsonI (j2007^ reported a mass of 2.7 Mq and classified this source as protostellar. 
iPhillips et all (|2QQlD classify this source (L BS18S) as pre-protostellar, however. 

— HOPS359: INutter fc Ward-ThompsonI (2007) classified this source as protostellar and measure a mass of 2.7 M©. 

— HOPS341: This source is strongly blended with HOPS340. We have listed some detections of the combined system 
in TableEl 

— HOPS354: lR.eipurth etld] (j2008D described this source in the context of the L1622 cloud. iBallv et al.l (j2009r ) also 
detected this source in their analysis of Spitzer IRAC images but did not analyze it in detail. 
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GALLERY OF IMAGES OF PAGS BRIGHT RED SOURGES 
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Fig. 15. — Same as Figure [3 showing 4^ x 4^ images of PBRs 061012. The IRAC-band emission associated with the source is clearly 
visible at 4.5 /im. Contours indicate the 870 /im emission levels at {0.1, 0.25, 0.4, 0.55, 0.7} Jybeam"-*^; the long wavelength sum-mm 
data trace the cold envelope material associated with the source. The 160 fim panel is shown with the lowest 870 //m emission contour. 




Fig. 16. — Same as Figure [3 showing 4' x 4' images and SED of PBRs 119019. Contours indicate the 870 /im emission levels at {0.1, 
0.2, 0.3, 0.4, 0.5} Jybeam"-*^. This source has prominent IRAC emission and is located in an IRAC 8 fim and MIPS 24 fim shadow. 
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Fig. 17. — Same as Figure showing 4' x 4^ images and SED of HOPS169. Contours indicate the 870 /xm emission levels at {0.1, 0.2, 
0.4, 0.6, 0.8} Jybeam-i. This source has clear outflow activity, traced by the IRAC emission, and appears to be at high inclination. 




Fig. 18. — Same as Figure [T] showing 4^ x 4^ images and SED of PBR source 019003. Contours indicate the 870 /im emission levels at 
{0.5, 1.0, 1.5, 2.0, 2.5} Jybeam"^. This source has indications outflow activity, as traced by the IRAC emission. The Herschel photometry 
may be strongly affected by blending due to the source location in a very dense fllament. 
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Fig. 19. — Top: Same as Figure[3 showing 5' x 5' images of three red sources: HOPS372 and 082012 (top) and 082005 (bottom). Contours 
indicate the 870 /xm emission levels at {0.25, 0.5, 0.75, 1.0, 1.25, 1.5} Jybeam-^. No IRAC emission is detected for 082005; however, this 
source is located in dense filamentary material traced by the sub-mm emission and an 8 /im absorption feature. 
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Fig. 20. Same as Figured showing 4' x 4' images and SED of PBRs 090003. Contours indicate the 870 /im emission levels at {0.25, 0.5, 
0.75, 1.0, 1.25, 1.5} Jybeam"-*^. The IRAC data show faint indi cations of extended em ission, possibly associated with an outflow activity 
or a second source. The 350 /im image and SED point are from IMiettinen et al.l (|2012l ) 
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Fig. 21. Same as FigureH 4' x 4^ images of PBRs 302002. Contours indicate the 870 /im emission levels at {0.25, 0.5, 0.75, 1.0, 1.25, 
1.5} Jy beam"-*^, tracing the cometary globule shaped region. The IRAC data show indications of emission associated with outflow activity; 
furthermore, the 4.5 /xm data show evidence that this source is observed at high inclination. This source is the second most massive source 
in our sample, with a best-flt M.env = 1.7 Mq. 
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Fig. 22. 

Jy beam"-*^. 



Same as Figure[7l 4' x 4' images of HOPS359. Contours indicate the 870 /im emission levels at {0.25, 0.5, 0.75, 1.0, 1.25, 1.5} 
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NGC2068 H0PS341 
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Fig. 23. — Same as Figure[3 4' x4' images of HOPS341. Contours indicate the 870 /im emission levels at {0.25, 0.5, 0.75, 1.0} Jybeam"^. 
The X -shaped morphology of this source in the IRAC bands indicated that tit is a binary. Indeed, the photometry of this source is strongly 
blended with HOPS340; nevertheless we include this source in our sample for completeness. 




Fig. 24. — Same as Figure [T] 4' x 4' images of PBRs 097002. Contours indicate the 870 /im emission levels at {0.25, 0.5} Jybeam"^. 




Fig. 25. — Same as Figure □ 4' x 4' images of HOPS354. Contours indicate the 870 /xm emission levels at {0.25, 0.5, 0.75, 1.0, 1.25} 
Jybeam-i. The IRAC images display clear indications of outflow activity; furthermore the source appears highly inclined. 
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TABLE 1 

Summary of HOPS Herschel PACS LI 641 observations 
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^ /I O OO o o 

— 7:49:29.88 


L1641 


1342218731 


;2) 


703 


8x8 


71 


5:40:40.150 


T. /I O . 1 O T1 

— 7:49:18.71 


L1641 


1342228163 


;4) 


844 


8x8 


320 


5:40:58.890 


T /I o OO o r 

— 7:48:02.05 


L1641 


1342228165 


;6) 


844 


5x5 


73 


5:40:42.910 


1-7 A e: oi oi 

— 7:45:01.91 


L1641 


1342228425 


;6) 


847 


5x5 


75 


r Ar\ c\ A /^orv 

5:40:24.620 


^7 A O OO O /'* 

— 7:43:08.26 


L1641 


1342227082 


;3) 


831 


5x5 


76 


5:40:26.090 


—7:37:32.02 


L1641 


1342205246 


;7) 


502 


5x5 


79 


5:39:57.200 


OO 1 O OO 

—7:30:19.89 


L1641 


1342205244 


;5) 


502 


8x8 


77 


5:40:44.670 


fr OO K A A 

—7:29:54.46 


L1641 


1342228427 


;8) 


847 


5x5 


123 


5:40:08.780 


—7:27:27.68 


L1641 


1342228161 


;2) 


844 


5x5 


26 


5:39:24.640 


1^ O/^ 1 O O 1 

— 7:26:13.81 


L1641 


1342218729 


;30) 


703 


8x8 


25 


5:39:55.200 


T. O /I . C O T1 

— 7:24:53.71 


L1641 


1342215589 


;9o) 


663 


8x8 


313 


5:39:33.300 


^ OO r" ^ o/^ 

— 7:22:57.36 


L1641 


1342227084 


;5) 


831 


5x5 


28 


5:38:56.470 


^ OO A A O O 

— 7:20:44.32 


L1641 


1342227086 


;7) 


831 


8x8 


30 


5:38:44.050 


^ "1 "1 yl O O O 

—7:11:49.89 


L1641 


1342204254 


;5) 


484 


8x8 


29 


5:39:06.710 


^ "1 "1 "1 o o o 

— 7:11:12.80 


L1641 


1342204252 


;3) 


484 


8x8 


32 


5:38:01.100 


r\ O^ 01 

—7:07:37.01 


L1641 


1342227045 


;6) 


830 


8x8 


312 


5:38:46.540 


— 7:05:37.46 


L1641 


1342205242 


;3) 


502 


5x5 


31 


r" o o /I /I o 

5:38:44.870 


^ OO o^ oo 

— 7:00:37.03 


L1641 


1342204256 


;7) 


484 


8x8 


33 


5:38:20.090 


r o o /I or 

—6:59:04.85 


L1641 


1342228171 


;2) 


844 


5x5 


35 


5:37:24.460 


r" O OO 1^1^ 

—6:58:32.77 


L1641 


1342227314 


;5) 


833 


5x5 


34 


5:37:59.990 


—6:57:27.50 


L1641 


1342205240 


;i) 


502 


8x8 


101 


5:37:17.090 


/I O /I o o o 

—6:49:49.33 


L1641 


1342227312 


;3) 


833 


5x5 


36 


5:37:52.390 


—6:47:18.67 


L1641 


1342227088 


;9) 


831 


8x8 


38 


5:36:22.050 


AV A O O 

—6:45:41.23 


L1641 


1342205238 


;9) 


502 


8x8 


40 


5:36:26.650 


oo o^ ^ A 

—6:38:27.74 


L1641 


1342227094 


;5) 


831 


8x8 


37 


5:37:17.280 


O/^ 1 O 1 o 

—6:36:18.18 


L1641 


1342227090 


;i) 


831 


8x8 


43 


V or' r' M MOM 

5:35:50.020 


r* ^ A r' o o 

—6:34:53.40 


L1641 


1342227310 


;i) 


833 


5x5 


50 


Oyl "ir' OOO 

5:34:15.880 


O yl o o ^o 

—6:34:32.70 


L1641 


1342217748 


;9) 


686 


8x8 


39 


5:36:41.330 


o /I oo oo 

—6:34:00.08 


L1641 


1342227092 


;3) 


831 


5x5 


41 


O/^ A A r\ 

5:36:19.440 


—6:29:06.79 


L1641 


1342227316 


;7) 


833 


5x5 


45 


5:35:34.120 


—6:26:41.70 


L1641 


1342215593 


;4) 


663 


8x8 


311 


J" O/I OM o/^r\ 

5:34:39.860 


or" 1/1 1 

—6:25:14.16 


L1641 


1342203649 


;5o) 


470 


5x5 


42 


J" O/^ OO /l/^O 

5:36:22.460 


r* o o o o 1 /f 

—6:23:39.14 


L1641 


1342205236 


;7) 


502 


8x8 


44 


5:36:36.980 


n -\ A r T oo 

—6:14:57.98 


L1641 


1342204258 


;9) 


484 


5x5 


51 


5:35:22.180 


1 O O/^ O /I 

—6:13:06.24 


L1641 


1342227318 


;9) 


833 


5x5 


47 


J" o/^ 1 ^ o/^r\ 

5:36:17.260 


—6:11:11.00 


L1641 


io4ZZZ i oZ4 




833 


5x5 


49 


5:35:52.000 


-6:10:01.85 


L1641 


1342227322 


;3) 


833 


5x5 


52 


5:35:33.210 


-6:06:09.65 


L1641 


1342227320 


;i) 


833 


5x5 


48 


5:36:31.360 


-6:01:16.81 


L1641 




;5) 


685 


8x8 


5 


5:35:07.960 


-5:56:56.40 


L1641 


1342204248 


;9) 


484 


8x8 


6 


5:35:24.560 


-5:55:33.42 


L1641 


1342227328 


;9) 


833 


5x5 


7 


5:36:19.020 


-5:55:25.46 


L1641 


1342227326 


;7) 


833 


5x5 


8 


5:35:04.400 


-5:51:00.76 


L1641 




;7) 


685 


8x8 


9 


5:33:30.710 


-5:50:41.03 


L1641 


1342217750 


;i) 


686 


8x8 


10 


5:36:10.100 


-5:50:08.34 


L1641 


1342227096 


;7) 


831 


5x5 


12 


5:34:46.830 


-5:42:28.72 


L1641 


1342204246 


;7) 


484 


8x8 


13 


5:35:17.340 


-5:42:14.51 


L1641 


1342227098 


;9) 


831 


5x5 


14 


5:34:30.440 


-5:37:47.44 


L1641 


1342204244 




484 


8x8 


^ Field center coordinates. 
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TABLE 2 

Summary of HOPS Herschel PACS ONC, NGC2024, NGC2068, and L1622 

OBSERVATIONS 



-— T 

HOPS 




R.A." 


; 


\ 

Region 


AOR ID 


— 

OD 


\ 

Map size 

X 


group name 


h:m:s 


o ./ .// 








15 


5:35:06.620 


—5:35:05.68 


ONC 


1342205234(5) 


502 


8x8 


308 


5:33:45.870 


-5:32:58.09 


ONC 


1342204433(4) 


487 


5x5 


16 


5:34:43.990 


—5:32:11.21 


ONC 


1342217448(9) 


685 


8x8 


17 


5:35:16.320 


—5:29:32.60 


ONC 


1342217450(1) 


685 


8x8 


18 


5:33:55.730 


—5:22:39.97 


ONC 


1342217752(3) 


686 


8x8 


200 


5:35:19.270 


—5:14:46.49 


ONC 


1342205232(3) 


502 


8x8 


130 


5:35:24.710 


— 5-09-06 02 


ONC 


1342205228(9) 


502 


8x8 


135 




— '^-nfi-S'^ 24 

O.VJU.OO.Zi^ 


ONC 




502 


8x8 


19 


5:35:23.300 


—5:00:35.73 


ONC 


1342204250(1) 


484 


8x8 


20 


5:35:13.770 


—4:54:57.33 


ONC 


1342217758(9) 


686 


8x8 


21 


5 •34- 32 340 


—4-53-54 26 


ONC 


1342217754f5') 


686 


8x8 


306 


5:35:32.280 


—4:46:48.47 


ONC 


1342191970(1) 


300 


5x5 


24 


5:35:23.340 


—4:40:10.45 


ONC 


1342217756(7) 


686 


8x8 


80 


5-40-51 710 


—2:26:48.62 


NGC2024 


1342226729(30) 


826 


5x5 


81 


5:41:28.940 


—2:23:19.36 


NGC2024 


1342226733(4) 


826 


5x5 


82 


5:41:23.740 


—2:16:51.10 


NGC2024 


1342228913(4) 


858 


8x8 


83 


5:41:42.180 


-2:16:26.20 


NGC2024 


1342227049(50) 


830 


8x8 


85 


5:42:02.620 


—2:07:45.70 


NGC2024 


1342226735(6) 


826 


5x5 


86 


5:41:43.560 


— 1:53:28.42 


NGC2024 


1342227047(8) 


830 


8x8 


89 


5:42:27.680 


— 1:20:01.00 


NGC2024 


1342205220(1) 


502 


5x5 


90 


5:43:04.370 


— 1:16:11.60 


NGC2024 


1342228376(7) 


849 


8x8 


91 


5:46:06.690 


—0:13:05.15 


NGC2068 


1342205218(9) 


502 


8x8 


92 


5:46:14.210 


-0:05:26.84 


NGC2068 


1342205216(7) 


502 


5x5 


93 


5:46:40.830 


+0:00:30.52 


NGC2068 


1342215363(4) 


662 


8x8 


94 


5:46:39.580 


+0:04:16.61 


NGC2068 


1342228365(6) 


848 


5x5 


302 


5:46:28.320 


+0:19:49.40 


NGC2068 


1342228374(5) 


849 


5x5 


303 


5:47:24.810 


+0:20:59.68 


NGC2068 


1342227966(7) 


843 


8x8 


96 


5:47:08.970 


+0:21:52.86 


NGC2068 


1342215587(8) 


663 


8x8 


128 


5:46:56.220 


+0:23:42.41 


NGC2068 


1342218727(8) 


703 


8x8 


301 


5:45:53.590 


+0:25:27.30 


NGC2068 


1342216450(1) 


675 


5x5 


97 


5:47:58.060 


+0:35:30.12 


NGC2068 


1342227969(70) 


843 


8x8 


98 


5:47:31.850 


+0:38:05.77 


NGC2068 


1342227971(2) 


843 


8x8 


300 


5:47:42.990 


+0:40:57.50 


NGC2068 


1342205214(5) 


502 


5x5 





5:54:15.240 


+1:43:15.59 


L1622 


1342215365(6) 


662 


8x8 


1 


5:54:55.370 


+1:45:03.08 


L1622 


1342218780(1) 


704 


8x8 


3 


5:54:23.540 


+1:49:17.78 


L1622 


1342218703(4) 


702 


8x8 


4 


5:54:36.260 


+1:53:54.00 


L1622 


1342218778(9) 


704 


8x8 



Note. — Note: the ONC field contains the extended Orion Nebula region, the NGC1977 
region, and OMC2/3; the NGC2024 field contains the NGC2024 HH region and the 
NGC2023 reflection nebula. The NGC2068 field includes the NGC2068 and NGC2071 re- 
flection nebulae as well as LBS23 region. 
^ Field center coordinates. 
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TABLE 3 

Herschel protostar candidate coordinates and photometry 



Source 


group 


R.A.« 


Decl.^ 


Field 


24 /im 


70 /im 


160 /im 


flag^ 


870 /im 


Tboi 


Lboi 




name 


h:m:s 


o ./.// 




[mjy] 


[mJy] 


[mJy] 




detection*^ 


[K] 


[L©] 


061012 


061 


05:42:48.87 


-08:16:10.70 


L1641 


<1.14 


703±35 


5634±845 


1 




32.1±0.9 


0.75±0.06 


119019 


119 


05:40:58.47 


-08:05:36.10 


L1641 


1.46±0.2 


1604±80. 


10745±1612 


1 


Yes 


34.4±0.8 


1.56±0.14 


026011 


026 


05:39:17.00 


-07:24:26.64 


L1641 


<1.30 


78.4±3.9 


2516±377 


1 


Yes 


23.0±1.0 


0.34±0.03 


313006 


313 


05:39:30.75 


—07:23:59.40 


L1641 


<1.25 


138±7. 


2198±330 


1 


Yes 


27.8±1.0 


0.28±0.03 


029003 


029 


05:39:13.15 


-07:13:11.69 


L1641 


1.19±0.2 


29.5±1.5 


365±55 


1 




42. Oil. 5 


0.06±0.01 


019003 


019 


05:35:23.92 


-05:07:53.46 


ONC 


<16.2 


2506±127 


31577±4737 


1 


Yes 


33.6±1.1 


3.16±0.31 


082005 


082 


05:41:29.40 


-02:21:17.06 


NGC2024 


<3.13 


506±25 


9308±140 


1 


Yes 


29.3±0.8 


1.02±0.11 


082012 


082 


05:41:24.94 


-02:18:08.54 


NGC2024 


6.51±0.4 


4571±229 


51254±7688 


1 


Yes 


32.2±0.9 


6.27±0.65 


090003 


090 


05:42:45.23 


-01:16:14.18 


NGC2024 


4.74±0.3 


3286±164 


16937±2541 


1 


Yes 


36.0±0.8 


2.71±0.24 


091015 


091 


05:46:07.65 


-00:12:20.73 


NGC2068 


<1.29 


648di32 


6615±992 


1 


Yes 


30.9±0.8 


0.81ib0.07 


091016 


091 


05:46:09.97 


-00:12:16.85 


NGC2068 


<1.14 


431di22 


5108±766 


1 


Yes 


29.1dz0.9 


0.65±0.06 


093005 


093 


05:46:27.75 


-00:00:53.81 


NGC2068 


<1.14 


1427ib71. 


12131ibl820 


1 


Yes 


30.8di0.9 


1.71ib0.15 


302002 


302 


05:46:28.24 


00:19:27.00 


NGC2068 


<1.14 


302dzl5 


7187dzl08 


1 


Yes 


28.6it0.9 


0.84dz0.09 


i~\ g-\ Pw i~\ i~\ «^ 

097002 


097 


05:48:07.76 


00:33:50.79 


NGG2068 


<1.14 


1049±52. 


7993±120 


1 


Yes 


33.4±0.9 


1.14±0.11 


300001 


300 


05:47:43.36 


00:38:22.43 


NGC2068 


<7.19 


478±24 


5042±756 


1 


Yes 


29.6±0.9 


0.65±0.06 


068006 


068 


05:41:11.79 


-07:53:35.09 


L1641 


9.19±0.2 


44.8±2.2 


146±22 


2 




149.7±5.06 


0.05±0.01 


038002 


038 


05:36:11.11 


-06:49:11.29 


L1641 


<1.25 


36.4±1.8 


537±81 


2 




40.4±1.0 


O.lOiO.Ol 


037003 


037 


05:37:00.35 


-06:37:10.95 


L1641 


8.85±0.2 


725±36 


3798±570 


2 




53.1±1.7 


0.43±0.03 


037008 


037 


05:37:34.31 


-06:35:20.33 


L1641 


<1.39 


27.7±1.4 


321±48 


2 




45.5±1.9 


0.04±0.01 


092011 


092 


05:46:26.17 


-00:04:45.31 


NGG2068 


2.31±0.2 


77.8±3.9 


363±54 


2 




104.1±5.63 


0.05±0.01 


093001 


093 


05:46:56.32 


-00:03:14.73 


NGC2068 


<1.14 


21.8±1.1 


217±33 


2 


Yes 


47.9±0.7 


0.05±0.01 


302004 


302 


05:46:16.55 


00:21:35.09 


NGC2068 


2.87±0.2 


108±5. 


283±42 


2 




54.8±1.8 


0.07±0.01 


096023 


096 


05:46:53.23 


00:22:10.05 


NGC2068 


<2.13 


147di7. 


1691di254 


2 




39.6dzl.4 


0.18ib0.02 


301003 


301 


05:46:02.15 


00:23:29.86 


NGC2068 


8.08ib0.2 


42.6ib2.1 


231di35 


2 




92.1di3.8 


O.OSdiO.Ol 


000011 


000 


05:54:32.10 


01:42:54.92 


L1622 


<1.25 


28.5±1.4 


147d=22 


2 




50.1di0.5 


0.06ib0.01 


053002 


053 


05:43:24.07 


-08:49:03.75 


kOri 


<1.38 


33.5ibl.7 


190±28 


3 




43.2ib0.3 


0.06ib0.01 


117004 


117 


05:41:40.40 


-08:41:40.60 


L1641 


1.17±0.2 


25.3±1.3 


234ib35 


3 




56.4±2.2 


0.04±0.01 


117014 


117 


05:41:29.28 


-08:36:14.60 


L1641 


1.20±0.2 


31.0±1.6 


237±36 


3 




32.2±2.4 


O.lliO.Ol 


119016 


119 


05:40:40.54 


-08:05:55.00 


L1641 


1.97±0.2 


42.7±2.1 


148±22 


3 




62.7±1.7 


0.05±0.01 


121011 


121 


05:41:37.90 


-07:55:44.35 


L1641 


1.46±0.2 


74.1±3.7 


193±29 


3 




51.6±1.5 


0.04±0.01 


025044 


025 


05:39:56.80 


-07:19:21.40 


L1641 


<1.18 


23.1±1.2 


150±22 


3 




42.2±0.2 


0.08±0.01 


030013 


030 


05:38:59.58 


-07:10:31.92 


L1641 


1.43±0.2 


25.9±1.3 


153±23 


3 




48. Oil. 


0.05±0.01 


031003 


031 


05:38:55.35 


-07:05:29.17 


L1641 


2.96±0.2 


49.9±2.5 


175±26 


3 




67.7±1.8 


0.06±0.01 


034010 


034 


05:37:54.76 


-06:56:59.65 


L1641 


<1.14 


61.3±3.1 


249±37 


3 




43.5±0.7 


0.15±0.01 


031037 


031 


05:38:28.20 


-06:56:40.16 


L1641 


<1.14 


64.3±3.2 


396±59 


3 




51.9±2.4 


0.06±0.01 


034014 


034 


05:38:14.92 


-06:53:03.56 


L1641 


1.89±0.2 


47.4±2.4 


235±35 


3 




65.5±3.5 


0.03±0.01 


036003 


036 


05:38:05.97 


-06:50:58.91 


L1641 


1.74±0.2 


67.4±3.4 


291±44 


3 




74.6±3.1 


0.04±0.01 


036006 


036 


05:37:45.07 


-06:50:02.39 


L1641 


2.83±0.2 


34.8±1.7 


143±22 


3 




60.1ibl.8 


0.04±0.01 


036011 


036 


05:37:42.72 


-06:47:08.31 


L1641 


2.68db0.2 


41.8di2.1 


159di24 


3 


Yes 


65.9dil.9 


O.OSdbO.Ol 


037011 


037 


05:37:37.80 


-06:34:43.35 


L1641 


2.42ib0.2 


58.0zb2.9 


409dz61 


3 




69.6zb3.0 


0.04ib0.01 


050006 


050 


05:34:23.01 


-06:32:58.00 


L1641 


2.82ib0.2 


45.5ib2.3 


313dz47 


3 




60.3dz2.3 


0.06±0.01 


037013 


037 


05:37:22.00 


-06:32:56.48 


L1641 


2.93±0.2 


52.0±2.6 


180±27 


3 




62.7±1.8 


0.05±0.01 


041001 


041 


05:36:28.68 


—06:30:42.13 


L1641 


2.72±0.2 


47.5±2.4 


262±39 


3 




90.7±4.8 


0.03±0.01 


009001 


009 


05:33:32.55 


-05:53:34.25 


ONC 


4.83±0.3 


74.9±3.8 


245±37 


3 




95.4±2.6 


0.05±0.01 


021010 


021 


05:34:19.63 


-04:53:23.54 


ONC 


11.3±0.5 


316±16 


678±10 


3 




63.8±2.2 


0.18±0.01 


303017 


303 


05:47:49.01 


00:20:26.47 


NGC2068 


4.64±0.2 


107±5. 


634±95 


3 




70.2±3.8 


O.lOiO.Ol 


303023 


303 


05:47:45.58 


00:21:14.68 


NGC2068 


4.36±0.2 


65.9±3.3 


234±35 


3 




104.1±4.81 


0.04±0.01 


098001 


098 


05:47:25.77 


00:33:37.43 


NGC2068 


1.73±0.2 


9.79±0.5 


231±35 


3 




43.5±2.1 


0.03±0.01 


097003 


097 


05:47:45.88 


00:34:12.71 


NGC2068 


1.37±0.2 


13.0±0.7 


137±21 


3 




78.1±4.4 


0.02±0.01 


000003 


000 


05:54:17.28 


01:40:18.68 


L1622 


8.63±0.4 


165±8. 


386±58 


3 




82.5±3.3 


O.lOiO.Ol 


000010 


000 


05:54:37.27 


01:42:52.39 


L1622 


5.10dz0.3 


64.7±3.2 


215±32 


3 




89.2dz3.9 


0.04dz0.01 


069001 


069 


05:40:46.20 


-08:04:35.12 


L1641 


5.13ib0.2 


325dzl6 


722±11 


4 




336.1ibl4.7 


0.29ib0.01 


026001 


026 


05:39:18.49 


-07:27:52.37 


L1641 


4.71di0.2 


223±11 


598±90 


4 




224.4dzll.2 


0.17±0.01 


306004 


306 


05:35:24.66 


-04:49:43.53 


ONC 


<124 


435ib22 


637ib96 


4 




119.1ib3.56 


0.33ib0.01 


006006 


006 


05:35:11.47 


-05:57:05.09 


L1641 


5.06±0.2 


126±6. 


4027±604 


5 


Yes 


27.7±0.6 


0.31±0.04 



Note. — Note: the ONC field contains the extended Orion Nebula region, the NGC1977 region, and OMC2/3; the NGC2024 field contains 
the NGC2024 HIT region and the NGC2023 reflection nebula. The NGC2068 field includes the NGC2068 and NGC2071 reflection nebulae as well 
as LBS23 region. 

^ PACS 70 /im source coordinates. 

^ This column indicates if a source is flagged as a reliable protostar (value = 1), if the source is considered less likely be a protostar (value = 2), 
if the source is flagged as extragalactic contamination (value = 3), if the SED shape is with IRAC photospheric emission (value = 4), and finally, 
one source has no IRAC coverage (value = 5). 
This column indicates if the source has a strong 870 /xm detection; sources with no data are either are assigned upper limits or have no coverage 
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TABLE 4 

Spitzer AND Herschel photometric properties of PACS bright red sources 



Source 


HOPS 


R.A.^ 


Decl.^ 


Field 


24 /im 


70 /im 


100 /im 


160 /im 


log 70/24^ 




group name 


h:m:s 


o ./.// 




[mJy] 


[mJy] 


[mJy] 


[mJy] 




061012 


061 


05:42:48.8 


-08:16:10.70 


L1641 


<1.14 


703±35 


2018±517 


5634±845 


>2.31 


119019 


119 


05:40:58.4 


-08:05:36.10 


L1641 


1.46±0.2 


1604±80. 


5789±148 


10745±1612 


2.56 


HOPS169 


040 


05:36:36.0 


-06:38:54.02 


L1641 


4.80±0.5 


5001±250 


15753±4033 


29975±4496 


2.54 


019003 


019 


05:35:23.9 


-05:07:53.46 


ONC 


<16.2 


2506±127 


4711±121 


31577±4737 


>1.71 


082005 


082 


05:41:29.4 


-02:21:17.06 


NGC2024 


<3.13 


506±25 


3003±769 


9308±140 


>1.73 


HOPS372 


082 


05:41:26.3 


-02:18:21.08 


NGC2024 


12. Oil. 2 


6178±309 


16217±4151 


31090±4664 


2.24 


082012 


082 


05:41:24.9 


—02:18:08.54 


NGC2024 


6.51±0.4 


4571±229 


20357±5211 


51254±7688 


2.37 


090003 


090 


05:42:45.2 


-01:16:14.18 


NGC2024 


4.74±0.3 


3286±164 


10914±2794 


16937±2541 


2.36 


HOPS358 


091 


05:46:07.2 


-00:13:30.86 


NGC2068 


422±42 


60681±3041 


104322±26706 


123207±18481 


1.68 


091015 


091 


05:46:07.6 


-00:12:20.73 


NGC2068 


<1.29 


648±32 


2543±651 


6615±992 


>2.22 


091016 


091 


05:46:09.9 


-00:12:16.85 


NGC2068 


<1.14 


431±22 


1977±506 


5108±766 


>2.10 


HOPS373 


093 


05:46:30.7 


-00:02:36.80 


NGC2068 


14.1±1.4 


5258±263 


20188±5168 


36724±5509 


2.10 


093005 


093 


05:46:27.7 


-00:00:53.81 


NGC2068 


<1.14 


1427±71. 


5373±138 


12131±1820 


>2.62 


302002 


302 


05:46:28.2 


00:19:27.00 


NGC2068 


<1.14 


302±15 


3101±794 


7187±108 


>1.95 


HOPS359 


303 


05:47:24.8 


00:20:58.24 


NGC2068 


22.8±2.3 


20758±1039 


43592±1116 


60409±9061 


2.48 


HOPS341 


128 


05:47:00.9 


00:26:20.76 


NGC2068 


14.0±1.4 


3001±301 


15138±3875 


25213±3782 


1.86 


097002 


097 


05:48:07.7 


00:33:50.79 


NGC2068 


<1.14 


1049±52. 


4163±107 


7993±120 


>2.49 


HOPS354 


000 


05:54:24.1 


01:44:20.15 


L1622 


28.9±2.9 


8492±426 


37423±9580 


39258±5889 


1.99 



Note. — Note: the ONC field contains the extended Orion Nebula region, the NGC1977 region, and OMC2/3; the NGC2024 field contains 
the NGC2024 HIT region and the NGC2023 refiection nebula. The NGC2068 field includes the NGC2068 and NGC2071 refiection nebulae 
as well as LBS23 region. 

^ Object coordinates are derived from the PACS 70 /im images, 
b log AFa70/AFa24 



TABLE 5 

APEX 350 AND 870 /iM photometry of PACS bright red sources 



Source 


350 /im 


350 /im" 


350 /im^ 


870 /im 


870 /im« 


870 /im' 




[Jy beam"-*^] 


[Jy] 


[Jy] 


[Jy beam"-*^] 


[Jy] 


[Jy] 


061012 


2.60 


3.93 


2.53 


<0.7 


<1.3 


<1.3 


119019 


3.38 


5.27 


3.63 


0.5 


0.9 


0.6 


HOPS169 








1.0 


1.6 


1.4 


019003 


8.29 


19.6 


7.18 


2.5 


5.0 


3.0 


082005 








0.8 


1.6 


1.0 


HOPS372 








<1.4 


<4.0 


<4.0 


082012 








2.7 


4.3 


3.6 


090003 


3.63^ 






1.7 


2.2 


1.9 


HOPS358 


13.0 


22.4 


17.3 


1.8 


3.3 


2.6 


091015 


2.14 


3.15 


2.26 


0.6 


1.3 


0.5 


091016 


2.50 


3.76 


2.81 


0.7 


1.2 


0.8 


HOPS373 


9.02 


12.9 


10.5 


1.5 


2.5 


2.0 


093005 


6.58 


10.1 


6.74 


1.4 


2.6 


1.6 


302002 








1.1 


1.7 


1.5 


HOPS359 








1.8 


2.9 


2.2 


HOPS341 








1.0 


1.5 


1.3 


097002 








0.6 


0.9 


0.8 


HOPS354 


13.8 


25.9 


18.7 


1.5 


2.7 


2.2 



^ Source fiux measured in an aperture with radius equal to 1 x FWHM, where the 
FWHM= 7.34'' and 19.0'' at 350 and 870 /im, respectively. 

^ Source fiux measured in the same aperture as (a) but with with local sky subtraction 
over radii equal to {1.5, 2.0} x FWHM. 
350 /im point from ^Miettinen et al.i (i2012|). 
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TABLE 6 

Orion PBRs observed properties: detections, environment, and previous detections 



vDKJ U-i V^t/ 


R.A. 


Decl. 


4 5 m 


Hip-h 
nign 


24 /im 


Note^ 






h.m.s 


o ./.// 


QGtGCtlOn 


incl " 


detection 




references 


061012 


05:42:48.8 


-08:16:10.7 


yes 




no 


ns 




119019 


05:40:58.4 


-08:05:36.1 


yes 




yes 






HOPS169 


05:36:36.0 


-06:38:54.0 


yes 


yes 


yes 


f,i 


Le88, Mo91, Za97, DaOO, St02, Jo06, Nu07, Ba09, Mel2 


019003 


05:35:24.0 


-05:07:50.1 


yes 




no 


f,c 


Me90, Ch97, Ni03, Ts03, Nu07, Mel2 


082005 


05:41:29.4 


-02:21:17.1 


no 




no 


f,i 


La91, La96, Mo99, Jo06, Nu07 


HOPS372 


05:41:26.3 


-02:18:21.1 


yes 




yes 


f,ns 


La91, La96, Mo99, Jo06, Nu07, Mel2 


082012 


05:41:24.9 


-02:18:08.5 


no 




yes 


f,ns 


La91, La96, Mo99, Wu04, Jo06, Nu07 


090003 


05:42:45.2 


-01:16:14.2 


yes 




yes 


i 


Mi09 


HOPS358 


05:46:07.2 


-00:13:30.9 


yes 




yes 


f,c 


St86, Li99, Re99, MiOl, Wu04, Nu07, Mel2 


091015 


05:46:07.7 


-00:12:19.1 


no 




no 


f,ns 


Li99, MiOl, Nu07 


091016 


05:46:10.0 


-00:12:15.4 


no 




no 


f,ns 


Li99, MiOl, Nu07 


HOPS373 


05:46:30.7 


-00:02:36.8 


yes 




yes 


f,i 


Ha83, La91, GiOO, MoOl, MiOl, Nu07, Mel2 


093005 


05:46:27.7 


-00:00:51.5 


yes 




no 


f,i 


La91, MoOl, MiOl, Nu07 


302002 


05:46:28.2 


+00:19:28.4 


yes 


yes 


no 


f,ns 


La91, JoOl, PhOl, MoOl, Sa03, Nu07 


HOPS359 


05:47:24.8 


+00:20:58.2 


no 




yes 




La91, MoOl, Nu07, Mel2 


HOPS341 


05:47:00.9 


+00:26:20.8 


yes 


yes 


yes 


b 


JoOl, MoOl, Nu07, SalO, Mel2 


097002 


05:48:07.7 


+00:33:52.7 


no 




no 


b/ns 




HOPS354 


05:54:24.2 


+01:44:20.1 


yes 


yes 


yes 


b 


Re08, Ba09, Mel2 



Note. — Le88: 'L evreaulti fl988); Mo91: 'M organ et al] ((19911); ZA97: Zavagno et al." (199?)- DaOO: 'Davis et al.J (120001 ): 
St02: Stanke et al. (2002^); Da02: Davis et al. (2009); Nu07: [Nutter & Ward-Thompson (2007) and (Di Francesco et al. 200 8]); Jo06: 
[Johnstone & Ballv (2006); Ts03: rTsuiimoto et al. (20031); Ch97: Chini et al. (1997); Ni03: Nielbock et al. (2003); Me90: Mezger et aP 
iTggiy): Mo 99: Mookeriea et al. (2009); La96: Launhardt et al. (1996); La91: Lada et al. (1991a); Wu04: Wu et al. (2004); Mi09: 
[ Miettinen et al. (2009); Re99 Reipurth et al. (1999); St86: Strom et al. (1986); Li99: Lis et al. (1999); MiOl: Mitchell et al. (2001); MoOl: 
[Motte et al. (2001); GiOO: Gibb & Little (2000); Ha83: Haschick et al. (1983); PhOl: Phillips et al. (2001); Sa03: Savva et al. (2001); 
JoOl: Johnstone et al. (2001); Re08: Reipurth et al. (2008); Ba09: Ballv et al. (2009); SalO: S adavov et al. (2010); Mel2: .Megeath et all 
(p012) 

^ Sources with indications of an observed high inclination orientation from the 4.5 fim images. 

^ ns = nearby source and indicates that we find a source close (< 30^^) to the target; f = filament and indicates that a significant 
level of elongation or extended emission is seen in the 870 /im data; i = isolated and indicates that no near-by source is observed; c = 
crowed, indicating that the source resides in region possibly with multiple sources, extended emission, or both and whose photometry 
may be affected by source blending; finally, b = binary, are sources that may have an unresolved secondary source and therefore blended 
photometry. 



TABLE 7 

Numbers and fractions of new sources and red PBRs found 
by region 



Region 


Field 


Total 


New 


New 


PBRs 


PBRs 






number 


number 


% 


number 


% 


All 




355 


56 


16% 


18 


5% 


Orion A 




274 


35 


13% 


4 


1% 




L1641^ 


200 


31 


16% 


3 


2% 




ONC 


74 


4 


5% 


1 


1% 


Orion B 




81 


21 


26% 


14 


17% 




NGC2024 


17 


3 


18% 


4 


24% 




NGC2068 


52 


15 


29% 


9 


17% 




L1622 


12 


3 


25% 


1 


8% 



^ kOti sources have been included here. 
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TABLE 8 

Orion PBRs and comparison Class sources: observed properties and modified black body fit 

parameters 



Source 


Tboi 


Lbol 


Lsmm /Li 1 

' DOl 




Lb 




peaK,iViiDiD 




[Kl 

L J 


L yJ i 


% 


[Kl 

L J 


L yJ i 


[M©] 




Ir^ J 


061012^ 


32.1±0.9 


0.75±0.06 


3.1 


19.2±0.4 


0.68±0.07 


0.20±0, 


,04 


133 


119019 


34.4±0.9 


1.56±0.14 


2.6 


19.9±0.4 


1.35±0.11 


0.32±0, 


.06 


127 


HOPS169 


35.4±0.9 


4.50±0.42 


2.4 


20.4±0.5 


3.85±0.41 


0.81±0, 


.20 


127 


019003^ 


33.6±1.1 


3.16±0.29 


1.1 


21.4±0.8 


1.52±0.30 


0.24±0. 


.18 


121 


082005^ 


29.3±0.8 


1.02±0.11 


5.4 


17.0±0.4 


1.02±0.13 


0.61±0. 


.15 


151 


HOPS372 


36.9±1.0 


4.90±0.45 


1.9 


21.3±0.6 


3.96±0.43 


0.65±0. 


.17 


121 


082012 


32.2±0.9 


6.27±0.65 


3.8 


18.6±0.4 


5.89±0.71 


2.20±0. 


.52 


139 


090003 


36.0±0.8 


2.71±0.24 


1.8 


21.4±0.4 


2.06±0.15 


0.33±0. 


.06 


121 


HOPS358 


44.3±0.9 


30.6±2.21 


0.6 


27.2±1.0 


19.5±1.23 


0.78±0. 


.14 


88. 


091015^ 


30.9±0.9 


0.81±0.07 


3.3 


18.7±0.4 


0.72±0.06 


0.24±0, 


.04 


136 


091016^ 


29.1±0.9 


0.65±0.06 


4.1 


17.9±0.3 


0.61±0.06 


0.26±0, 


.05 


142 




Qd nj_n n 


K onj_n A Q 
O.ZUzLU.4o 


O Q 
Z.O 


on oj_n A 
zU.zzLU.4 


A Ol J_n QO 


0.93±0, 


.16 


1 0'7 

iz ( 


093005^ 


30.8±0.9 


1.71±0.15 


3.4 


18.8±0.4 


1.57±0.15 


0.52±0. 


.10 


136 


302002^ 


28.6±0.9 


0.84±0.10 


6.1 


16.2±0.3 


0.82±0.11 


0.63±0, 


.15 


157 


HOPS359 


39.3±0.9 


12.6±1.00 


1.3 


23.4±0.7 


9.56±0.84 


0.91±0. 


.23 


115 


HOPS341 


36.3±1.1 


3.62±0.37 


3.1 


19.1±0.5 


3.16±0.35 


0.97±0. 


.24 


133 


097002^ 


33.4±0.9 


1.14±0.11 


2.8 


19.4±0.5 


0.99±0.11 


0.27±0. 


.07 


130 


HOPS354 


37.4±0.8 


7.53±0.74 


1.9 


21.0±0.4 


5.84±0.45 


1.05±0. 


.19 


124 


B335^ 


40.3±1.3 


0.55±0.04 


2.6 


20.0±0.7 


0.47±0.04 


O.lliO, 


,02 


130 


CB 130^ 


55.2±3.3 


0.22±0.01 


10.9 


12.7±0.5 


0.17±0.01 


0.55±0, 


,11 


199 


CB 244-SMMl^ 


63.2±2.6 


1.72±0.14 


1.8 


20.3±0.8 


l.lOiO.ll 


0.24±0, 


,04 


127 


CB68^ 


49.2±1.4 


0.68±0.06 


1.9 


19.8±1.1 


0.42±0.05 


O.lOiO, 


,02 


130 


VLA1623-243-^ 


34.0±1.2 


2.82±0.32 


2.4 


21.3±0.5 


2.55±0.24 


0.48±0, 


,09 


124 



^ Our Lsmm /Lbol may be underestimated as the Lsmm values are measured from the modified black-body 
fits to the SEDs (see text). 

^ Best-fit modified black-body parameters; the parameters shown here are derived from fitting the A > 
70 fim SED points. Excluding the 70 /im point and fitting only the A > 100 /im points systematically 
increases the derived masses by 20% on average, decreases the temperature by 5% on average, and decreases 
the luminosity by 7% on average. 

^ The ma sses derived here will increase by a factor of if we assume an ISM-type dust model 
(|Draine Lee 1984). Furthermore, we assume a gas-to-dust ratio of 110 (see text). 

^ The values of T^qI and Lbol for these sources should be considered upper limits because we include the 

24 fim upper limits in our analysis for these sources. 

^ Additional well-studied isolated Class sources from ILaunha rdt et al.' ("2012*) and IStutz et al.l (|2Q10I ) 

are shown for comparison; the distances to these sources lie between 100 pc and 250 pc. 

^ SED from J. Green and DIGIT team, private communication, 2012, and Green et al., in prep. 



